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1. Introduction 

1.1 Historical overview 

In 1930 the Austrian physicist Wolfgang Pauli proposed the existence of particles called 
neutrinos, denoted as i^, as a "desperate remedy" to account for the missing energy in a 
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type of radioactivity called beta decay. At the time physicists were puzzled because nuclear 
beta decay appeared to violate energy conservation. In beta decay, a neutron in an unstable 
nucleus transforms into a proton and emits an electron, where the radiated electron was 
found to have a continuous energy spectrum. This came as a great surprise to many 
physicists because other types of radioactivity involved gamma rays and alpha particles 
with discrete energies. Pauli deduced that some of the energy must have been taken away 
by a new particle emitted in the decay process, the neutrino, which carries energy and has 
spin 1/2, but which is massless, electrically neutral and very weakly interacting. Because 
neutrinos interact so weaJily with matter, Pauli bet a case of champagne that nobody would 
ever detect one, and they became known as "ghost particles". Indeed it was not until a 
quarter of a century later, in 1956, that Paiili lost his bet and neutrinos were discovered 
when Clyde Cowan and Fred Rcincs detected antineutrinos emitted from a nuclear reactor 
at Savannah River in South Carolina, USA. 

Since then, after decades of painstaking experimental and theoretical work, neutri- 
nos have become enshrined as an essential part of the accepted quantum description of 
fundamental particles and forces, the Standard Model (SM) of particle physics. This is 
a highly successful theory in which elementary building blocks of matter are divided into 
three generations of two kinds of particle - quarks and Icptons. It also includes three of the 
fundamental forces of Nature, the strong (g), electromagnetic (7) and weak {W,Z) forces 
carried by spin 1 force carrying bosons (shown in parentheses) but does not include gravity. 
There are six flavours of quarks. The leptons consist of the charged electron e~ , muon fi" 
and tau r~, together with three electrically neutral particles - the electron neutrino Vg, 
muon neutrino and tau neutrino Vr which are our main concern here. 

The first clues that neutrinos have mass came from an experiment deep underground, 
carried out by an American scientist Raymond Davis Jr., detecting solar neutrinos [1]. It 
revealed only about one-third of the number predicted by theories of how the Sun works 
pioneered by John Bahcall [1]. The result puzzled both solar and neutrino physicists. 
However, some Russian researchers, Mikheyev and Smirnov, developing ideas proposed 
previously by Wolfenstein in the U.S., suggested that the solar neutrinos might be chang- 
ing into something else. Only electron neutrinos are emitted by the Sun and they could 
be converting into muon and tau neutrinos which were not being detected by the experi- 
ments. This effect called "neutrino oscillations" , as the types of neutrino interconvert over 
time from one kind to another, was first proposed some time earlier by Pontecorvo [2]. 
The precise mechanism for "solar neutrino oscillations" proposed by Mikheyev, Smirnov 
and Wolfenstein involved the resonant enhancement of neutrino oscillations due to matter 
effects, and is known as the MSW effect [3]. 

Neutrino oscillations are analogous to coupled pendulums, where oscillations in one 
pendulum induce oscillations in another pendulum. The coupling strength is defined in 
terms of something called the "lepton mixing matrix" U which relates the basic Standard 
Model neutrino states, Vg, V/^, Ut-, associated with the electron, muon and tau, to the 
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neutrino mass states vi, 1^2, and 1^3 with (real and positive) masses mi, m2, and ms [4], 
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According to quantum mechanics it is not necessary that the Standard Model states z^e, i^fi, 
Vt be identified in a one-one way with the mass eigenstates z^i, f2, and 1/3, and the matrix 
elements of U give the quantum amplitude that a particular Standard Model state contains 
an admixture of a particular mass eigenstate. The probability that a particular neutrino 
mass state contains a particular SM state may be represented by colours as in Fig. |. Note 
that neutrino oscillations are only sensitive to the differences between the squares of the 
neutrino masses Am?^- = mf — m|, and gives no information about the absolute value of 
the neutrino mass squared eigenvalues mf. There are basically two patterns of neutrino 
mass squared orderings consistent with the atmospheric and solar data as shown in Fig. |l]. 
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Figure 1: The probability that a particular neutrino mass state contains a particular SM state 
may be represented by colours as shown in the key. Note that neutrino oscillation experiments 
only determine the difference between the squared values of the masses. Also, while rn^ > it is 
presently unknown whether m| is heavier or lighter than the other two, corresponding to the left 
and right panels of the figure, referred to as normal or inverted mass squared ordering, respectively. 
Finally the value of the lightest neutrino mass (sometimes referred to as the neutrino mass scale) 
is presently unknown and is represented by a question mark in each case. 
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As with all quantum amplitudes, the matrix elements of U are expected to be complex 
numbers in general. The lepton mixing matrix U is also frequently referred to as the 
Maki-Nakagawa-Sakata (MNS) matrix C/mns [4], and sometimes the name of Pontecorvo 
is added at the beginning to give C/pMNS- The standard parameterisation of the PMNS 
matrix in terms of three angles and at least one complex phase, as recommended by the 
Particle Data Group (PDG) [5], will be discussed later. 

Before getting into details, here is a quick executive summary of the implications of 
neutrino mass and mixing following from Fig. |l]: 

• Lepton flavour is not conserved, so the individual lepton numbers Lg, L^, are 
separately broken 

• Neutrinos have tiny masses which are not very hierarchical 

• Neutrinos mix strongly unlike quarks 

• The SM parameter count is increased by at least 7 new parameters (3 neutrino masses, 
3 mixing angles and at least one complex phase) 

• It is the first (and so far only) new physics beyond the SM 

The idea of neutrino oscillations was first confirmed in 1998 by the Japanese experi- 
ment Super-Kamiokande (SK) [6] which showed that there was a deficit of muon neutrinos 
reaching Earth when cosmic rays strike the upper atmosphere, the so called "atmospheric 
neutrinos". Since most neutrinos pass through the Earth unhindered, Super-Kamiokande 
was able to detect muon neutrinos coming from above and below, and found that while 
the correct number of muon neutrinos came from above, only about a half of the expected 
number came from below. The results were interpreted as half the muon neutrinos from 
below oscillating into tau neutrinos over an oscillation length L of the diameter of the 
Earth, with the muon neutrinos from above having a negligible oscillation length, and so 
not having time to oscillate, yielding the expected number of muon neutrinos from above. 

In 2002, the Sudbury Neutrino Observatory (SNO) in Canada spectacularly confirmed 
the "solar neutrino oscillations" [7]. The experiment measured both the flux of the elec- 
tron neutrinos and the total flux of all three types of neutrinos. The SNO data revealed 
that physicists' theories of the Sun were correct after all, and the solar neutrinos Ve were 
produced at the standard rate but were oscillating into z/^ and z^t, with only about a third 
of the original Ue flux arriving at the Earth. 

Since then, neutrino oscillations consistent with solar neutrino observations have been 
seen using man made neutrinos from nuclear reactors at KamLAND in Japan [8] , and neu- 
trino oscillations consistent with atmospheric neutrino observations have been seen using 
neutrino beams fired over hundreds of kilometres as in the K2K experiment in Japan [9], 
the Fermilab-MINOS experiment in the U.S. [10] or the CERN-OPERA experiment in Eu- 
rope. Further long-baseline neutrino beam experiments are in the pipeline, and neutrino 
oscillation physics is entering the precision era, with superbeams and a neutrino factory on 
the horizon. 
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Following these results several research groups showed that the electron neutrino has 
a mixing matrix element of | Ue2 1 ~ 1/ which is the quantum amplitude for to contain 
an admixture of the mass eigenstate 1^2 corresponding to a massive neutrino of mass 1712 ~ 
0.008 electronvolts (eV) or greater (where -y/ rrig — mf ~ 0.008 eV). By comparison the 
electron has a mass of about half a megaelectronvolt (MeV). Put another way, the mass 
state 1^2 contains roughly equal probabilities of I'e, and v^- sometimes called trimaximal 
mixing, corresponding to the three equal red, green and blue colours associated with 
in Fig. |. The muon and tau neutrinos were observed to contain approximately equal 
amplitudes of the third neutrino v-^ of mass , | 1 w | Urz I ~ 1/ where a normalised 
amplitude of l/\/2 corresponds to a 1/2 fraction of 1^3 in each of i/^ and Vr, leading to a 
maximal mixing and oscillation of ■v^ Ur- Put another way, the mass state v-^ contains 
roughly equal probabilities of and Uj- called maximal mixing, corresponding to the two 
equal green and blue colours associated with m\ in Fig. ||. Interestingly, the value of ms 
is not determined and it could be anywhere between zero and 0.3 eV, depending on the 
mass scale and ordering. Although at least one neutrino mass must be 0.05 eV or greater 
(where y^Ittt-I — mW ~ 0.05 eV), this could be either or m2, as shown in Fig. [l|. 

According to the early results from the CHOOZ nuclear reactor experiment [11], the 
electron neutrino v^, could only contain a very small amount of the third neutrino mass 
eigenstate ^3, \Ue-i\ < 0.2. Evidence for non-zero C/gS was first provided by T2K, MINOS 
and Double Chooz [12]. Recently the Daya Bay [13], RENO [14], and Double Chooz [15] 
collaborations have measured \Ue2,\ ~ 0.15. Put another way, the mass state 1^3 has a 
probability of containing V(, of about (0.15)^, corresponding to the small amount of red 
colour associated with m| in Fig. 1. As we shall see, this element being non-zero excludes 
a number of simple mixing patterns and models which were previously proposed, and has 
led to a number of new developments. 

1.2 Where we stand 

The main experimental milestones from 1998-2012 may be summarised as follows: 

• 1998 - SK confirms that atmospheric are converted to another neutrino type, 
probably consistent with near maximal mixing \U^'i\ ~ [f^rsl ~ l/v^ 

• 2002 - SK, SNO and the older neutrino experiments such as Homestake and the 
Gallium experiments results are combined in a global fit pointing towards the large 
(but non-maximal) mixing and conversion of solar neutrinos in the core of the Sun 

• 2002 - SNO confirms that solar Uf, are converted to a linear combination of and Ur 
with approximate trimaximal mixing \Ue2\ ^ \U^2\ ^ |f^T2| ^ l/v^ 

• 2004 - Reactor antineutrinos z7e are observed by KamLAND to oscillate with a prob- 
ability consistent with the solar neutrino oscillations 

• 2006 - Accelerator neutrinos from Fermilab are observed over a long baseline 
(LBL) at MINOS with a disappearance probability consistent with the atmospheric 
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oscillation results, providing a high precision confirmation of a similar observation 
from KEK to SK (K2K) in 2004 

• 2010 - LBL accelerator neutrinos from CERN appear at OPERA as Vr 

• 2011 - T2K and MINOS observe and excess of accelerator neutrinos appearing as 
z/e, consistent with non-zero C/e3 

• 2012 - Daya Bay, RENO and Double Chooz observe the disappearance of reactor 
antineutrinos and measure jt/esl ~ 0.15 

1.3 The challenges ahead for experiment 

Despite the above observations, neutrinos remain the least understood particles. Of the (at 
least) 7 new parameters which must be present due to neutrino mass and mixing, only 5 
are currently measured, namely the three mixing angles and two mass squared differences. 
For example none of the CP violating phases are currently measured, although there are 
plans to measure one of these phases in next generation neutrino oscillation experiments. 
However, since the neutrino oscillations are only sensitive to mass squared differences, the 
lightest neutrino mass cannot be measured by oscillation experiments. Also the present 
experiments are not sensitive enough to uniquely determine the ordering of the neutrino 
square masses m\, m\, m|, although it is known that m\ > The neutrino mass scale 
(i.e. the mass of the lightest neutrino) is not known, although, as discussed later, there are 
cosmological reasons to believe that none of the neutrino masses can exceed about 0.3 eV. 
Hence the lightest neutrino mass should be somewhere between zero and 0.3 eV. However 
cosmology is not sensitive to whether neutrino mass is of the Dirac or Majorana kind.^ 
In principle, if neutrinoless double beta decay were observed, it could simultaneously be 
used to measure both the lightest neutrino mass and show that it is Majorana (and future 
measurements could shed light on the additional phases associated with Majorana masses). 
However, despite our ignorance, we know that neutrino masses are much smaller than the 
other charged fermion masses, and this already represents something of a puzzle. 

From the experimental perspective, the main known unknowns of neutrino mass and 
mixing may be summarised as: 

• The neutrino mass squared ordering (normal or inverted) 

• The neutrino mass scale (i.e. the mass of the lightest neutrino, presumably between 
zero and 0.3 eV) 

• The nature of neutrino mass (Dirac or Majorana) 

• The CP violating phase measurable in neutrino oscillations (the so called Dirac 
phase 5, although it is also present if neutrino mass is Majorana) 

• The two possible further CP violating phases associated with Majorana neutrino 
masses (not present if neutrino mass is Dirac) 

^See Subsection |l.5| for the definition of Dirac and Majorana neutrino masses. 
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Neutrino physics has now entered the precision era, at least as far as the measured param- 
eters are concerned. T2K is presently running [16] and will provide accurate measurements 
of the atmospheric neutrino mass squared difference and mixing angle, while NOz/A [17], 
presently under construction, will provide complementary information about the mass or- 
dering. Future neutrino oscillation experiments, under discussion [18], will give more ac- 
curate information about the mass squared splittings Am^^- = mf — rnj, mixing angles, the 
mass squared ordering (commonly but incorrectly referred to as the "mass hierarchy"), and 
the neutrino mass scale (i.e. the mass of the lightest neutrino mass eigenstate, which will 
indeed decide if neutrino masses involve a significant mass hierarchy). The ultimate goal 
of oscillation experiments, however, is to measure the so far undetermined CP violating 
oscillation phase 5 and there is considerable activity in this area [19] to determine the best 
way to do this. 

1.4 The nature and scale of neutrino mass 

Oscillation experiments are not by themselves capable of telling us anything about the 
nature or mass scale of neutrino mass. They can, however, shed light on the neutrino mass 
ordering as mentioned above. There are basically four ways to elucidate the mysteries 
surrounding neutrino masses. 

1. Neutrinoless double beta decay experiments (for a recent review see e.g. [20]) effec- 
tively measure the 1-1 element of the Majorana neutrino mass matrix corresponding 



and can validate the Majorana nature of neutrinos. There was a claim of a signal 
in neutrinoless double beta decay corresponding to m^^ = 0.11 — 0.56 eV at 95% 
C.L. [21]. However this claim was criticised by two groups [22], and in turn this 
criticism has been refuted [23]. Experiments such as GERDA should report soon and 
decide this question [20]. 

2. Oscillation experiments can measure the sign of Am^2 resolve normal from in- 
verted mass squared orderings. 

3. Independently of whether neutrinos are Dirac or Majorana, the Tritium beta decay 

experiment KATRIN [24] will tell us about the absolute scale of neutrino mass down 
to about 0.35 eV. Such experiments measure the "electron neutrino mass" defined by 



4. More model dependently, cosmology can in principle probe the sum of neutrino 
masses, and hence the lightest neutrino mass mughtest) down to very small values [25]. 
In future detection of energetic neutrinos from gamma ray bursts, neutrino telescopes 
could also provide important astrophysical information, and may provide another 
means of probing neutrino mass, and even quantum gravity [26]. 
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Figure 2: vs. miightost: The red and light red regions represent the model independent values 
that the inverted neutrino mass ordering can take based on the central value and la deviation of 
a recent global fit of neutrino parameters. The blue and light blue regions are the analogue of this 
for the normal neutrino mass ordering. The gold regions correspond to the golden ratio prediction 
for 171013 in both the normal and inverted orderings, resulting from the ^5 inverse mass sum rule. 



In Fig. ^ we show the allowed range of the effective mass parameter for neutrinoless 
double beta decay, mfsp, see Eq. (|1.2|), as a function of the lightest neutrino mass mughtest 
for both the normal (blue) and inverted (red) neutrino mass squared orderings consistent 
with the one sigma range of parameters taken from a recent global fit as discussed in [27] 
from which this figure is taken. ^ Also shown (gold) is the restricted region in a model 
based on a discrete family symmetry which involves the golden ratio (GR) [27]. Such 
restricted regions follow from relations between the neutrino masses which can generally 
arise in models based on discrete family symmetry. For example, the neutrino masses may 
be related by a sum rule of the form, 

ami + /3m2 = m^, (1-4) 

where a, /3 are model dependent constants. If the model involves a seesaw mechanism, then 
the right-handed neutrino masses may be similarly related leading to inverse relationships 
between light physical neutrino masses of the form, 

^ + ^ = ^. (1.5) 

Till 1712 

where 7, 6 are model dependent constants.^ In certain models analogous relations may 
arise with the neutrino masses replaced by their square roots. All these mass sum rules 

^Fig. ^ is generated using a code whose original version was developed in [28]. 

^Note that 5 here is nothing to do with the CP violating phase denoted by the same Greek letter. 
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have been recently studied in [29]. The Ar, GR model [27] mentioned above involves an 
inverse sum rule as in Eq. (|1.5| ) with 7 = (5 = 1. This may be compared to the A4 model 
in [30] which involves an inverse mass sum rule with 7 = 1 and 6 = —2, or the A(96) 
model in [31] with 7 = 1 and 5 = ib2i, and so on. The appearance of the complex constant 
reminds us that in all the (inverse) mass sum rules there are CP violating phases associated 
with Majorana neutrino masses which are implicit. 

1.5 The origin of neutrino mass 

It is worth recalling why the observation of non-zero neutrino mass and mixing is evidence 
for new physics beyond the SM. The most intuitive way to understand why neutrino mass 
is forbidden in the Standard Model, is to understand that the Standard Model predicts 
that neutrinos always have a "left-handed" spin - rather like rifle bullets which spin counter 
clockwise to the direction of travel. In fact this property was first experimentally measured 
in 1958, two years after the neutrino was discovered, by Maurice Goldhaber, Lee Grodzins 
and Andrew Sunyar. More accurately, the "handedness" of a particle describes the direction 
of its spin vector along the direction of motion, and the neutrino being "left-handed" means 
that its spin vector always points in the opposite direction to its momentum vector. The 
fact that the neutrino is left-handed, written as vl, implies that it must be massless. 
If the neutrino has mass then, according to special relativity, it can never travel at the 
speed of light. In principle, a fast moving observer could therefore overtake the spinning 
massive neutrino and would see it moving in the opposite direction. To the observer, the 
massive neutrino would therefore appear right-handed. Since the Standard Model predicts 
that neutrinos must be strictly left-handed, it follows that neutrinos are massless in the 
Standard Model. It also follows that the discovery of neutrino mass implies new physics 
beyond the SM, with profound implications for particle physics and cosmology. 
Neutrinos are massless in the Standard Model for three independent reasons: 

• There are no right-handed neutrinos i^/j 

• There are only Higgs doublets (and no Higgs triplets) of SU (2) l 

• There are only renormalisable term 

In the SM, the three massless neutrinos fg, 1^^, Vt are distinguished by separate lepton 
numbers Lg, L^, Lr- Neutrinos and antineutrinos are distinguished by total conserved 
lepton number L = + + Lr- To generate neutrino mass we must relax one or more 
of the above three conditions. For example, by adding right-handed neutrinos the Higgs 
mechanism of the Standard Model can give neutrinos the same type of mass as the Dirac 
electron mass or other charged lepton and quark masses, which would generally break the 
separate lepton numbers L^, L^, L^, but preserve the total lepton number L. However it is 
also possible for neutrinos to have a new type of mass of a type first proposed by Majorana, 
which would also break L. 

From the theoretical perspective, the main unanswered question is the origin of neu- 
trino mass, and in particular the smallness of neutrino mass. The simplest possibility is 
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that neutrinos have Dirac mass just hke the electron mass in the SM, namely due to a 
term hke yoLHiyji, where L is a lepton doublet containing vl, H is a Higgs doublet and 
ur is a right-handed neutrino. The observed smallness of neutrino masses implies that the 
Dirac Yukawa coupling yo must be of order 10~^^ to achieve a Dirac neutrino mass of 
about 0.1 eV. Advocates of Dirac masses point out that the electron mass already requires 
a Yukawa coupling of about 10~^, so wc are used to such small Yukawa couplings. In 
this case, all that is required is to add right-handed neutrinos lyji to the SM and we are 
done. Well, almost. It still needs to be explained why the i^r have zero Majorana mass, 
after all they are gauge singlets and so nothing prevents them acquiring (large) Majorana 
mass terms Mrrvrvji where Mrr could be as large as the Planck scale. Moreover, Majo- 
rana masses offer a unique (and testable) way to generate neutrino masses (since neutrinos 
do not carry electric charge) even without right-handed neutrinos. The simplest way to 
generate Majorana mass is via yM^LL where A is a Higgs triplet and yM is a Yukawa 
coupling associated with Majorana mass. Alternatively, at the effective level, Majorana 
neutrino mass can result from some additional dimension 5 operators which couple two 
lepton doublets L to two Higgs doublets H first proposed by Weinberg [32] , 

—HL^kHL, (1.6) 

where k has dimension [mass]~^. This is a non-renormalisable operator, so it violates 
one of the tenets of the SM. In order to account for a neutrino mass of order 0.1 eV 
requires k ~ 10^-*^^ GcV^^. This suggests a new high energy mass scale M in physics, a 
small dimensionless coupling associated with k, or both. There are basically five different 
proposals for the origin of neutrino mass: 

• The seesaw mechanisms [33-35] (Weinberg operator e.g. from large Majorana mass 
M = Mjiji for right-handed neutrinos i/r) 

• i2-parity violating supersymmetry [36] (Weinberg operator from TeV scale Majorana 
mass for neutralinos %) 

• TeV scale loop mechanisms [37, 38] (Majorana mass from extra Higgs doublets and 
singlets at the TeV scale) 

• Extra dimensions [39] (Dirac mass with small yo due to right-handed neutrinos ur 
in the bulk) 

• String theory [40,41] (new mechanisms for generating large Majorana mass for right- 
handed neutrinos i^r from Planck or string scale physics) 

These different mechanisms axe reviewed in [42]. In this review we shall mainly focus on 
the seesaw mechanism which may be incorporated into a theory of flavour. 

It has been one of the long standing goals of theories of particle physics beyond the 
Standard Model to predict quark and lepton masses and mixings. With the discovery of 
neutrino mass and mixing, this quest has received a massive impetus. Indeed, perhaps 
the greatest advance in particle physics over the past decade has been the discovery of 
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neutrino mass and mixing involving large mixing. The largeness of the lepton mixing angles 
contrasts with the smallness of the quark mixing angles, and this observation, together with 
the smallness of neutrino masses, provides new and tantalising clues in the search for the 
origin of quark and lepton flavour. For example, it is amusing to note that the smallest 
lepton mixing may be related to the largest quark mixing, \Ue3\ ~ 6c /\/^ where 9c is the 
Cabibbo angle. The quest to understand the origin of the three families of quarks and 
leptons and their pattern of masses and mixing parameters is called the flavour puzzle, 
and motivates the introduction of family symmetry. In particular, as we shall see, lepton 
mixing provides a motivation for discrete family symmetry, which will form the central 
part of this review. As we shall also see, such theories demand a high precision knowledge 
of the lepton mixing angles, beyond that currently achieved. 

1.6 About this review 

It should be mentioned at the outset that there are good and fairly up to date reviews 
already in the literature, for example: [43-45], although only the last one was written 
after Daya Bay and RENO. It should be remarked that the signal of another independent 
mass splitting from the LSND accelerator experiment [46] would either require a further 
light neutrino state with no weak interactions (a so-called "light sterile neutrino" ) or some 
other non-standard physics. This effect has not been confirmed by a similar experiment 
KARMEN [47], and a subsequent experiment MiniBooNE [48] has not decisively resolved 
the issue. Since there is no solid evidence for light sterile neutrinos, in this review we shall 
not discuss this subject any further, but refer the interested reader to a recent dedicated 
discussion in [49]. Instead, in this review, we shall exclusively focus on the three active 
neutrino paradigm. 

The starting point for the present review is the one that was written by one of us about 
ten years ago [50]. At that time it had just become apparent, after the first SNO results 
in 2002, that the solar mixing angle was large, which together with the large atmospheric 
mixing angle, meant that there were two large mixing angles in the lepton sector. The solar 
mixing being large effectively killed many neutrino mass models that had previously been 
proposed consistent with small solar mixing. This was actually the second great extinction 
of models, the first being after the discovery of a large atmospheric mixing angle by Super- 
Kamiokande in 1998. A decade after the last review, in 2012 we are now in an analogous 
position, namely that Daya Bay and RENO have just measured the reactor angle and 
shown it to be quite sizeable, killing neutrino mass models consistent with a very small 
reactor angle. Just as the review article in [50] was written shortly after the second great 
extinction in 2002, so the present review article is being written just after the third great 
extinction in 2012, so once again it is an opportune moment to identify new and surviving 
model species which may come to dominate the theory landscape over the coming years. 
In particular the fate of discrete family symmetry, which was to some extent motivated 
by the possibility of the reactor mixing being zero, or very small, will be fully discussed. 
We emphasise that we are now in the unique position in the history of neutrino physics of 
knowing not only that neutrino mass is real, and hence the Standard Model at least in its 
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minimal formulation is incomplete, but also we finally have information on all three mixing 
angles. 

As in the previous review, we focus on theoretical approaches to understanding neu- 
trino masses and mixings in the framework of the seesaw mechanism, assuming three active 
neutrinos. The goal of such models is to account for two very large mixing angles, and 
one Cabibbo-sized mixing angle, and a pattern of neutrino masses consistent with obser- 
vation. We give a strong emphasis to classes of models where large mixing angles can 
arise naturally and consistently with a neutrino mass hierarchy. We show that if one of 
the right-handed neutrinos contributes dominantly in the type I seesaw mechanism to the 
heaviest neutrino mass, and a second right-handed neutrino contributes dominantly to the 
second heaviest neutrino mass, then large atmospheric and solar mixing angles may be 
interpreted as simple ratios of Yukawa couplings. This is of course the sequential domi- 
nance (SD) mechanism [51-53]. Sequential dominance is not a model, it is a mechanism 
in search of a model. The conditions for sequential dominance, such as ratios of Yukawa 
couplings being of order unity for large mixing angles, and the required pattern of right- 
handed neutrino masses are put in by hand and require further theoretical input such as 
family symmetry. It is interesting to note that, without further constraints, SD generically 
has long predicted hierarchical neutrinos with a normal ordering and large reactor angle of 
order |f/e3| ~ 0{m2/m^) ~ 0.2. However, in order to achieve more precise predictions, SD 
needs to be combined with family symmetry. 

The use of discrete family symmetry was mainly motivated by the hypothesis of exact 
tri-bimaximal (TB) mixing [54] defined by: 



The TB mixing and discrete family symmetry approach gained much impetus over the 
past decade. Given that the measurement of the reactor mixing [C/esj ~ 0.15 by Daya 
Bay and RENO kills exact TB mixing, the obvious question is what is the impact on the 
discrete family symmetry approach, which was largely inspired by TB mixing. This timely 
question will be addressed by the present review. There are a huge number of proposals in 
the literature, but not so many surviving the measurement of the reactor angle. The simple 
answer to the question is the discrete family symmetry approach is alive and kicking after 
Daya Bay and RENO. However, theorists have been forced to work harder to go beyond 
the simple mixing pattern of TB mixing which is now excluded. The simple discrete family 
symmetries proposed to account for TB mixing may still be viable as a leading order (LO) 
approximation, but higher order (HO) corrections may play a more important role than 
anticipated in many models. Alternatively, perhaps larger finite groups are relevant where 
the LO approximation already predicts a non-zero reactor angle. Yet another possibility 
is that the discrete family symmetry may be implemented indirectly, as in the sequential 
dominance approach, in new ways. All these interesting possibilities will be discussed. 



\Ue2,\ 
|f^M3| 
\Ue2\ 



0, 

\Urz\ 
\U,.2\ 



\Ur2\ = 1/V3. 



(1.7) 
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The layout of the remainder of the review article is as follows. In Section |2| we in- 
troduce and review the current status of neutrino masses and mixing angles. We also 
parameterise the PMNS mixing matrix in two different ways, whose equivalence is dis- 
cussed in Appendix In Section ^ we discuss patterns of lepton mixing that have been 
proposed, starting with simple mixing patterns such as bimaximal (BM), tri-bimaximal 
(TB), bi-trimaximal (BT) and golden ratio (GR) mixing. These all may apply to the 
neutrino mixing, which is then corrected by charged lepton mixing corrections, to give 
acceptable PMNS mixing, leading to solar mixing angle sum rules. The closeness of the 
TB mixing pattern to the data suggests a parametrisation of the PMNS matrix in terms of 
deviations from TB mixing, which is also introduced. Using these deviations, several TB 
variants are introduced and discussed, including tri-bimaximal-reactor (TBR) mixing and 
trimaximal (TM) mixing in two forms, namely where the first and second columns of the 
PMNS matrix take the TB values, called TMl and TM2, respectively. Section ^ is devoted 
to the seesaw mechanisms, which are central to this review, in both the simplest versions, 
called the type I, and including other types II and III, as well as alternative versions. We 
show how the type I seesaw mechanism may be applied to the hierarchical case in a very 
natural way using sequential dominance. Section |5| contains a mini-review of finite group 
theory and may be skipped by those readers who are already familiar with this subject. In 
Section |6| we give a pedagogical introduction to discrete family symmetry, and its direct or 
indirect implementation in model building. Section ^ is devoted to the direct model build- 
ing approach in which different subgroups of the discrete family symmetry are preserved 
in the neutrino and charged lepton sectors, and discusses the associated vacuum align- 
ments arising from the breaking of the discrete family symmetry using flavons. We also 
discuss the model building strategies following Daya Bay and RENO. Section ^ contains 
an analogous discussion for the indirect approach in which the discrete family symmetry is 
completely broken by flavons, but special vacuum alignments lead to particular mixing pat- 
terns, including new viable patterns with a non-zero reactor angle. In Section ^ we briefly 
review grand unifled theories (GUTs) such as SU{5) and how they may be combined with 
discrete family symmetry in order to account for all quark and lepton masses and mixing. 
In Section 10 we discuss three model examples which combine an SU{5) GUT with the 
discrete family symmetries A4, 6*4 and A(96). Section 11 concludes the review. We also 
present appendices dealing with more technical issues which may provide useful model 
building tools. Appendix ^ proves the equivalence between different parametrisations of 
the neutrino mixing matrix and gives a useful dictionary. Appendix |^ gives the full three 
family neutrino oscillation formula in terms of deviations from TB mixing. Appendix ^ 
catalogues the generators and Clebsch-Gordan coefficients of A4, 84, and Tj. 



2. Neutrino masses and mixing angles 



The history of neutrino oscillations dates back to the work of Pontecorvo who in 1957 [2] 
proposed — )• P oscillations in analogy with K ^ K oscillations, described as the mixing of 
two Majorana neutrinos. Pontecorvo was the first to realise that what we call the "electron 
neutrino" for example is really a linear combination of mass eigenstate neutrinos, and that 
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this feature could lead to neutrino oscillations of the kind ^ v^. Later on MSW proposed 
that such neutrino oscillations could be resonantly enhanced in the Sun [3]. The present 
section introduces the basic formalism of neutrino masses and mixing angles, and gives an 
up-to-date summary of the current experimental status of this fast moving field. 

2.1 Three neutrino mixing ignoring phases 

The minimal neutrino sector required to account for the atmospheric and solar neutrino 
oscillation data consists of three light physical neutrinos with left-handed flavour eigen- 
states, Ue, v^, and defined to be those states that share the same electroweak doublet 
as the corresponding left-handed charged lepton mass eigenstates. Within the framework 
of three- neutrino oscillations, the neutrino flavour eigenstates fg, f^, and Vr are related to 
the neutrino mass eigenstates vi, U2, and ^3 with mass mi, m2, and m^, respectively, by a 



3x3 unitary matrix called the lepton mixing matrix C/pMNS introduced in Eq. (LI). 

Assuming the light neutrinos are Majorana, C/pMNS can be parameterised in terms of 
three mixing angles 9ij and three complex phases Sij . A unitary matrix has six phases but 
three of them are removed by the phase symmetry of the charged lepton masses. Since 
the neutrino masses are Majorana there is no additional phase symmetry associated with 
them, unlike the case of quark mixing where a further two phases may be removed. 

To begin with, let us suppose that the phases are zero. Then the lepton mixing matrix 
may be parametrised by a product of three Euler rotations, 

t^PMNS = -R23-R13-R12, (2.1) 




R 



13 



^ Ci3 Si3 

1 I , Ri2= -si2 C12 I , (2.2) 

Y -Sl3 Ci3 




with Cij = COS 9ij and Sij = sin 6ij . Note that the allowed range of the angles is < 9ij < ^ . 

Ignoring phases, the relation between the neutrino flavour eigenstates z/g, f^, and I'r 
and the neutrino mass eigenstates i^i, 1^2, and z/3 is therefore given as a product of three 
Euler rotations in Eq. ( |2.1| ) as depicted in Fig. |^. 

2.2 Atmospheric neutrino mixing 

In 1998, the Super-Kamiokande experiment published a paper [6] which represents a wa- 
tershed in the history of neutrino physics. Super-Kamiokande measured the number of 
electron and muon neutrinos that arrive at the Earth's surface as a result of cosmic ray 
interactions in the upper atmosphere, which are referred to as "atmospheric neutrinos". 
While the number and angular distribution of electron neutrinos is as expected, Super- 
Kamiokande showed that the number of muon neutrinos is significantly smaller than ex- 
pected and that the flux of muon neutrinos exhibits a strong dependence on the zenith 
angle. These observations gave compelling evidence that muon neutrinos undergo flavour 
oscillations and this in turn implies that at least one neutrino has a non-zero mass. The 
standard interpretation is that muon neutrinos are oscillating into tau neutrinos. 
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Figure 3: The relation between the neutrino flavour eigenstates v^., v^^., and Vr and the neutrino 
mass eigenstates vi, 1/2, and 1/3 in terms of the three mixing angles O12, O13, 623- 

As a first approximation, one can set the reactor angle ^13 to zero, and assume that 
I ^"^-32 1 ^ I ^"^21 1 • Current atmospheric neutrino oscillation data can then approximately 
be described by simple two-state mixing, 

_ ( cos 6*23 sin 6*23 ) ( ^2 

Vr ) \- sin 6*23 cos 6*23 j 

and the two-state oscillation formula describing the probability that a converts to a i^t, 

P(i/^ ^ Vt) = sin^ 26*23 srn^ {1.27 /^ml^L/E), (2.4) 




where 

Amfj = mj — rrij , (2.5) 

and nii are the physical neutrino mass eigenvalues associated with the mass eigenstates z/j. 
Am|2 is in units of eV^, the baseline L is in km and the beam energy E is in GeV. 

The data can be well described by approximately maximal atmospheric mixing | C/^3 1 ^ 
I Ut-3 I ~ 1/ \/2 ■ This corresponds to 



sm c/23 



I/V2, (2.6) 



which means that the angle ^23 ~ 7r/4 or 45° and we identify the heavy "atmospheric 
neutrino" of mass ma as being approximately 

.3-^^. (2.7) 
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2.3 Solar neutrino mixing 



Super-Kamiokandc was also sensitive to the electron neutrinos arriving from the Sun, the 
"solar neutrinos" , and independently confirmed the reported deficit of such solar neutrinos 
long reported by other experiments. For example Davis's Homcstake Chlorine experiment 
which began data taking in 1970 consisted of 615 tons of tetrachloroethylene, and uses 
radiochemical techniques to determine the ^^Ar production rate. The SAGE and Gallex 
experiments contained large amounts of ^^Ga which is converted to ^^Ge by low energy elec- 
tron neutrinos arising from the dominant pp reaction in the Sun. The combined data from 
these and other experiments implied an energy dependent suppression of solar neutrinos 
which was interpreted as due to flavour oscillations. Taken together with the atmospheric 
data, this required that a second neutrino has a non-zero mass. The standard interpreta- 
tion is that the electron neutrinos disappear due to an oscillation formula which involves 
a "solar neutrino" of mass m2 given approximately by 



consistent with trimaximal solar mixing \Ue2\ ~ \Ufj,2\ ~ \Ur2\ ~ l/\/3- This corresponds 
to 



or 012 ~ 35°. 

SNO measurements of charged current (CC) reaction on deuterium were sensitive ex- 
clusively to while the neutral current (NC) reaction as well as the elastic scattering (ES) 
off electrons were also sensitive to and I'r- The neutrino flux derived from the CC reac- 
tions was significantly smaller than the one obtained from NC and ES. This immediately 
disfavoured oscillations of f e to sterile neutrinos which would lead to a diminished flux of 
electron neutrinos, but equal CC, NC and ES fluxes. On the other hand the observations 
were consistent with oscillations of i^e to active neutrinos ly^ and i^r since this would lead to 
a larger NC and ES rate. The SNO analysis was nicely consistent with both the hypoth- 
esis that electron neutrinos from the Sun oscillate into other active flavours, and with the 
Standard Solar Model prediction. The results from SNO including the data taken with salt 
inserted into the detector to boost the efficiency of detecting the NC events [7], strongly 
favoured the large mixing angle (LMA) MSW solution. In other words, after SNO, there 
was no longer any solar neutrino problem: we had instead solar neutrino mass m2^- 

KamLAND was a more powerful reactor experiment that measured z/g produced by 
surrounding nuclear reactors. KamLAND observed a signal of neutrino oscillations over 
the LMA MSW mass range, confirming the LMA MSW region "in the laboratory" [8]. 
KamLAND and SNO results when combined with other solar neutrino data especially that 
of Super-Kamiokande uniquely specify the LMA MSW [3] solar solution with three active 
light neutrino states, and approximately trimaximal solar mixing. KamLAND has thus not 
only confirmed solar neutrino oscillations, but has also uniquely specified the LMA solar 
solution, heralding a new era in neutrino physics. 




(2.9) 
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2.4 Reactor neutrino mixing 

Until recently, the reactor angle ^13 was not measured, only limited by CHOOZ, a reactor 
experiment that failed to see any signal of neutrino oscillations over the Super-Kamiokande 
mass range. CHOOZ data from Vf. — )• Vf. disappearance not being observed provided a 
significant constraint on ^13 over the Super-Kamiokande preferred range of Am|2 [H]: 

sin^ 26I13 < 0.2. (2.10) 

Direct evidence for ^13 was first provided by T2K, MINOS and Double Chooz [12]. Re- 
cently the Daya Bay [13], RENO [14], and Double Chooz [15] collaborations have measured 
sin2(2^i3): 

Daya Bay : sin2(26'i3) = 0.089 ± 0.010(stat.) ± 0.005(syst.) , 
RENO : sin2(26'i3) = 0.113 ± 0.013(stat.) ± 0.019(syst.) , (2.11) 
Double Chooz : sin2(26'i3) = 0.109 ± 0.030(stat.) ± 0.025(syst.) . 

This corresponds to 

|i7e3| =sin% «0.15, (2.12) 

or a reactor angle ^13 ~ 9°. 



2.5 Three neutrino mixing including phases 

If the reactor angle were zero then there would be no CP violation in neutrino oscillations. 
The measurement of the reactor angle means that we cannot ignore the presence of phases 
any more. Including the phases, assuming the light neutrinos are Maj orana, t^pMNS can 
be parameterised in terms of three mixing angles ^jj, a Dirac phase 5, together with two 
Maj orana phases /3i,/32, as follows [5], 



C^PMNS = R23U13R12P12, 



(2.13) 



where 



U 



13 



Cl3 




sise-'^' 

1 



\ -si3e'^ ci3 



Pi 



12 




(2.14) 



and R23 and R12 were defined below Eq. (2.1), giving, 



C^PMNS 



C12C13 
-S12C23 - Cl2Sl3S23e' 



S12C13 



5136 



C12C23 - Sl2'Sl3'S23e'" C13S23 I -Pl2- (2.15) 
\ S12S23 - Cl2Sl3C23e'^ -C12S23 - •Sl2Sl3C23e*^ C13C23 



Alternatively the lepton mixing matrix may be expressed as a product of three complex 
Euler rotations [55], 

t^PMNS = U23U13U12, (2-16) 
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parameter 


Forero et al 


Fogli et al 


Gonzalez-Garcia et al 


Am|i[10"^eV2] 


7.62 ± 0.19 


'•^^-0.22 


7.50 ±0.185 


Amii[10-3eV2] 


n + 06 


9 40+0.06 

^■^•J-o.io 

„ .„+o 11 


9 47+0.069 
^•^'-0.067 
„ -o+O 042 


sin^ 6^12 


0.320l°;g}? 


0.307l°:°l^ 


0.30 ±0.013 


sin^ 6^23 


6no+°-°26 

U.DUU_Q 


0.386l°:°^t 


41+0-037 

'-'•^-'^-0.025 
n 41+0.037 0^ r, c;q+0.021 


sin^ 6113 


U.UZ4D_o.o028 
U.UZOU_Q QQ27 


0.0241 ± 0.0025 
0.0244l°:»i 


0.023 ± 0.0023 


6 


(0.80 ± l)7r 
-(0.03± l)7r 


(1.08l°J«)vr 
(1.09i°:i)7r 


(1.67l°;^^)7r 



Table 1: Neutrino oscillation parameters summary. For Am§]^, sin^023, sin^^ia, and 6 the upper 
(lower) row corresponds to normal (inverted) neutrino mass ordering. The best fit values and la- 
errors are shown. The subtleties associated with these numbers are discussed in the respective 
references Forero et al [56], Fogli et al [57] and Gonzalez-Garcia et al [58]. In particular [58] quotes 
two different global fits, depending on the assumptions made about reactor fluxes, where we only 
quote the first ("free flux"). Furthermore, the precise definition of the atmospheric neutrino mass 
splitting Am|]^ differs slightly between the three global fits. 
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(2.17) 



(2.18) 



(2.19) 



The equivalence of different parametrisations of the lepton mixing matrix, and the 
relation between them is discussed in [52] with the results summarised in Appendix If 
the neutrinos are Dirac, then the phases /3i = /32 = 0, but the phase 6 remains. 



2.6 Global fits 

In Table |l| we give the global fits of the neutrino mixing parameters. For Am|^, sin^ ^23; 
sin^ ^13 , and 5 the upper (lower) row corresponds to normal (inverted) neutrino mass 
ordering. The best fit values and la errors are shown from Forero et al [56] , Fogli et al [57] 
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Figure 4: The mixing angles obtained from the three global fits [56-58]. The upper three panels 
correspond to the results for normal neutrino mass ordering (NO), while the lower three panels are 
for an inverted mass ordering (10). Shown are the best fit values (green) as well as the la (red) 
and 3cr (blue) intervals. Note that the solar angle is insensitive to the mass ordering. 



and Gonzalez-Garcia et al [58]. The results for the mixing angles are graphically contrasted 
in Fig. |[ 

For the normal mass ordering, we shall take the average values and errors which ap- 
proximately encompass the one sigma ranges of all three global fits (ignoring the solution 
of 023 in the second octant found by Forero et al [56]): 



sin2 6'i2 = 0.31 ±0.02, 
sin^ 6*23 = 0.41 ± 0.05, 



sm &is 



0.024 ± 0.003. 



(2.20) 
(2.21) 
(2.22) 



These values may be compared to the tri-bimaximal predictions sin^ 6 



12 



0.33, sin^ 6*23 



0.5 and sin^ ^13 = 0, showing that TB mixing is excluded by the reactor angle. Alternatively 
we may write, remembering that these are one sigma ranges in the squares of the sines and 
not the sines themselves, 

sin 012 = 0.56 ±0.02, (2.23) 

sin 6'23 = 0.64 ±0.04, (2.24) 

sin6li3 = 0.155 ±0.01. (2.25) 

In terms of the angles themselves we have, approximately, in round figures, 

012 = 34° ± 1°, (2.26) 

023 = 40° ± 3°, (2.27) 

6*13 = 9° ±0.5°. (2.28) 
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A few comments are relevant about these angles. Firstly the errors are not linear, since, 
for one thing, the global fits are made in terms of the squares of the sines of the angles. 
Having said this, there is a preference for the atmospheric angle to be in the first octant 
(i.e. less than 45°) and hence not maximal mixing. Secondly, as already noted, the solar 
angle is still consistent with trimaximal mixing (i.e. 35.26°) but there is a preference for it 
to be slightly smaller. 



3. Patterns of lepton mixing and sum rules 
3.1 Simple forms of neutrino mixing 

Below we give three examples of simple patterns of mixing in the neutrino sector which all 
have si3 = and S23 = C23 = l/\/2- Inserting these values in Eq. ( p.lD we obtain a PMNS 
matrix of the form. 



( C12 S12 

S12 C12 1 
V2 

\ £12. _£12 J_ 
\ \/2 \/2 \/2 , 



(3.1) 



where the zero subscript reminds us that this form has ^13 = (and ^23 = 45°). 

For golden ratio (GR) mixing [59], the solar angle is given by tan0i2 = l/</), where 
(j) = {\ + \/5)/2 is the golden ratio which implies 612 = 31.7°. There is an alternative 
version where cos ^12 = 0/2 and 612 = 36° [60], which we refer to as GR' mixing. 

For bimaximal (BM) mixing (see e.g. [61-63] and references therein), WG insert 512 — 
ci2 = l/\/2 (012 = 45°) into Eq. (O), 




For tri-bimaximal (TB) mixing [54], alternatively we use S12 
?12 = 35.26°) in Eq. (|3l|), 



(3.2) 



1/^/3, 



C12 



C/tb 



V 



1 1 1 

yM V2 



(3.3) 



These simple examples of neutrino mixing are now all excluded by the data. However 
they may still play a role in model building and we will revisit them when we consider 
charged lepton and other corrections below. 



3.2 Deviations from tri-bimaximal mixing 

From a theoretical or model building point of view, one significance of the global fits is that 
they exclude the tri-bimaximal lepton mixing pattern [54] in which the solar mixing angle 
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is trimaximal, the atmospheric angle is maximal and the reactor angle vanishes. When 
comparing global fits to TB mixing it is convenient to express the solar, atmospheric and 
reactor angles in terms of deviation parameters (s, a and r) from TB mixing [64]: 

1 1 r 

sin 6*12 = -^(1 + s), sin 6*23 = -^(1 + o), ^ ^'^'^'^ 

A related expansion is given in [65]. From these definitions we can write, 

s = \/3sin0i2 — 1, a = \/2sin023 — 1, r = \/2sin0i3. (3-5) 

Using this last form, and Eqs. (|^,|]2|j2^), we find the following values and ranges for 
the TB deviation parameters: 

s = -0.03 ±0.03, a = -0.10 ±0.05, r = 0.22 ± 0.01, (3.6) 

assuming a normal neutrino mass ordering. As well as showing that TB is excluded by the 
reactor angle being non-zero, Eq. (3.6) shows a preference (at the two sigma level) for the 



atmospheric angle to be below its maximal value and also a slight preference (at the one 
sigma level) for the solar angle to be below its trimaximal value. In general, this parametri- 
sation shows that the solar angle must be quite close to trimaximal, while the atmospheric 
angle may be far from bimaximal, with the reactor angle necessarily very far from zero. In 
any expansion in terms of these parameters, it should be a good approximation to work to 
first order in s, and possibly a, although it is worth working to second order to obtain the 
most accurate results. In Appendix |B| the PMNS matrix is expanded to second order in 
r,s,a, and the neutrino oscillation formulae including matter effects are given to a similar 
level of approximation (results taken from [64] ) . Note that the global fit values in Eq. ( |3.6| ) 
are consistent with, 

s = -5A^ a = -A\/2, r = R\, (3.7) 

where A is the Wolfenstein parameter and S, A, R are numbers all of order unity. In fact 
present data is consistent with S = A = R = l8X the one sigma level. In the next 
subsection we consider the simple case where s = a = and r = A which is within the two 
sigma range. 

3.3 Tri-bimaximal-Cabibbo mixing 

The recent data is consistent with the remarkable relationship [66], 

sin 61(7 A 

where A = 0.2253 ± 0.0007 [5] is the Wolfenstein parameter. The above ansatz implies a 
reactor angle of 

^i3«-^«9.2°, (3.9) 

where Oq ~ 13° is the Cabibbo angle. One may combine the relation in Eq. ( |3.8D with TB 
mixing to yield tri-bimaximal-Cabibbo (TBC) mixing [67]: 

.13 = ^, .12 = ^, .23 = ^. (3.10) 
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This corresponds to s = a = and r = A and leads to the following approximate form of 
the mixing matrix [67], 



C^TBC 



^'i2 + 0(A3), (3.11) 



corresponding to the mixing angles, 

^13 « 9.2°, ^12 = 35.26°, 023 = 45°. (3.12) 

This is consistent with the data at the two sigma level since it ignores the preference for 
the atmospheric angle to be in the first octant. 

3.4 Charged lepton contributions to neutrino masses and mixing angles 

The mixing matrix in the lepton sector, the PMNS matrix ?7pMNS) is defined as the matrix 
which appears in the electroweak coupling to the W bosons expressed in terms of lepton 
mass eigenstates. With the mass matrices of charged leptons Mg and neutrinos m^^ written 

4 

as 

C = -eZMeSR - ^iTZmlLi^l + H.c. , (3.13) 
and performing the transformation from flavour to mass basis by 

Me Vj^ = diag(me, m^, m^), "i^^ Vj^ = diag(mi, m2, ms), (3.14) 
the PMNS matrix is given by 

C/PMNS = Ve^V^^ . (3.15) 

Here it is assumed implicitly that unphysical phases are removed by field redefinitions, and 
C^PMNS contains one Dirac phase and two Majorana phases. The latter are physical only 
in the case of Majorana neutrinos, for Dirac neutrinos the two Majorana phases can be 
absorbed as well. 

As shown in [52,68-71] the lepton mixing matrix can be expanded in terms of neutrino 
and charged lepton mixing angles and phases to leading order in the charged lepton mixing 
angles which are assumed to be small, 

5236-*^^^ « s-^^e-''^-^ - eh4,e-'''^^ , (3.16) 
0136-'^'^ ^ ei^e-'^'^^ - 6'f3C^3e-^^i3 - e^2S23e'^"^2''3"^^2^ , (3.17) 
s^.2e-'''^ ^ sl^e-''"^^ + ei^s'i^e-^''^^ + ei^c'i^s'i^e'^'"^^-''^^'^ - ei^d^^c'i^e-^'r. , (3.I8) 

where we have dropped the subscripts L for simplicity. Clearly ^13 receives important 
contributions not just from O^^, but also from the charged lepton angles 0f2, and 6\^. In 



* Although we have chosen to write a Majorana mass matrix, all relations in the following are independent 
of the Dirac or Majorana nature of neutrino masses. 
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models where is extremely small, 613 may originate almost entirely from the charged 
lepton sector. Charged lepton contributions could also be important in models where 
^12 ~ since charged lepton mixing angles may allow consistency with the LMA MSW 
solution. 

Note that it is useful and possible to be able to diagonalise the mass matrices ana- 
lytically, at least to first order in the small 13 mixing angle, but allowing the 23 and 12 
angles to be large, while retaining all the phases. The procedure for doing this is discussed 
for a hierarchical and an inverted hierarchical neutrino mass matrix in [52]. The analytic 
results enable the separate mixing angles and phases associated with each of the unitary 
transformations Ve^ and vjj^ to be obtained in many useful cases of interest. 



3.5 Solar mixing sum rules 

In many models the neutrino mixing matrix has a simple form Uq in Eq 
S23 = C23 = l/\/2 and s'is 



|3J| ), where 

0, while the charged lepton mixing matrix has a CKM-like 



structure, in the sense that is dominated by a 12 mixing, i.e. its elements (14^)13, 
(Ki)23, (^5^)31 and (14^)32 are very small compared to (14^) 12 and {Ve^)2i, where in 
practice we take them to be zero. In this case we are led to a solar sum rule [69-71] derived 
from C/pMNS = Kj^f^Oi which takes the form, 
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(3.19) 



The important point to notice is that the 3-1, 3-2 and 3-3 elements of f/pMNS in Eq. ( 3.19| ) 
are uncorrected by charged lepton corrections and are the same as those of C/q, and also the 
1-3 element of f/pMNS has a simple form. By comparing Eq. ( 3.19|) to the PDG parametri- 
sation of ?7pmns in Eq. (p. 15 ) we find the relations. 
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(3.20) 
(3.21) 
(3.22) 
(3.23) 
(3.24) 



Using Eq. we see that, to leading order in ^13, the atmospheric angle is unchanged 

from its maximal value by the assumed form of the charged lepton corrections. To this 
approximation, it is then straightforward to expand these results to obtain the more useful 
approximate form of the sum rule [69-71], 
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+ 6*13 COS 5. 



(3.25) 
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Given the accurate determination of the reactor angle in Eq. ( |2.28| ) (^13 ~ 9° ± 0.5°) and 
the solar angle Eq. (|23^ ) (6*12 ^ 34° ± 1°) the sum rule in Eq. ( |3.25| ) yields a favoured 
range of cos 5 for each of the cases e'(2 = 35.26°, 45°, 31.7°, 36° for the cases of TB, BM, 
GR, GR', namely cos 5 — 0.2, —1, 0.2, —0.2, or cos 5 ^ — A,— 1,A,— A, respectively. For 
example, for TB neutrino mixing, the sum rule in Eq. ( pS] ) may be written compactly as. 



s ^ r cos 5. 



(3.26) 



-A^ from r ~ A, we need to have 



In order to obtain the values in Eq. (3/?), namely s 
cos 5 ~ —A. 

This approach relies on a Cabibbo-sized charged lepton mixing angle as is clear from 
Eq. ( 3.20 ) which, together with Eq. ( p.8[ ), shows that we need ~ A in order to account 
for the observed reactor angle, starting from one of the simple classic patterns of neutrino 
mixing. This is not straightforward to achieve in realistic models [67,72], which would 
typically prefer smaller charged lepton mixing angles such as ~ A/3. This suggests 
that the neutrino mixing angle is not zero, but has some non-zero value closer to the 
observed reactor angle. In the next subsection we consider this possibility. 



3.6 Atmospheric mixing sum rules 

In this subsection we consider simple alternative patterns related to TB mixing which allow 
a non-zero reactor angle. When looking for variants of TB mixing, it is useful to start from 
the general expansion around TB mixing in Eq. (|3.4D [64], which to leading order, gives a 



PMNS mixing matrix, as in Eq. (B.l), 
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Clearly TB mixing in Eq. ( |3.3D corresponds to s = a = r = 0. If we set s = a = but 
retain a non-zero value of r then this defines tri-bimaximal-reactor mixing (TBR) [73] , 
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(3.28) 



This is very constrained, in particular it requires maximal atmospheric mixing a = 0. We 
can maintain trimaximal (TM) mixing defined by s = but relax maximal atmospheric 
mixing, allowing both a non-zero a and r. 
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There are two interesting special cases of TM mixing in which the first or second column of 
the mixing matrix reduce to those of the first or second column of the TB mixing matrix, 
referred to as TMl and TM 2 mixing, namely, 
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and 



a = r cos 5, 
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(3.30) 



(3.31) 



(3.32) 



(3.33) 



a = — (r/2) cos 5. 

These TMl and TM2 relations, both with s = 0, are examples of atmospheric sum rules to 
first order in A. In order to obtain the values in Eq. (^), namely a ^ —A/2 with r w A, 
we see that TMl predicts cos 5 ~ — 1/2 and TM2 predicts cos 6 ~ 1. 

The above atmospheric sum rules are valid to leading order in A. The exact TMl 
relations (for both a and s) are obtained by equating PMNS elements to the first column 
of the TB mixing matrix: 
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(3.36) 



where these lead to Eq. ( 3.31| ) when expanded to leading order. 

The exact sum rule relations for TM2 are obtained by equating PMNS elements to the 
second column of the TB mixing matrix: 



S12C13 
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(3.37) 
(3.38) 
(3.39) 



where these lead to Eq. ( 3.33| ) when expanded to leading order. 

From Eqs. ( ^ , ^!37|) we see that to leading order in S13 the solar angle is unchanged 
from its TB value for both TMl and TM2, corresponding to s = as discussed earlier, but 
to second order in S13, the solar angle deviates and this deviation is different for TMl and 
TM2. 
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4. The seesaw mechanisms 

The starting point of the seesaw mechanisms is the dimension 5 operator in Eq. (|1.6|) which 
we repeat below, 

-^HL^kHL. (4.1) 

One might wonder if it is possible to simply write down an operator by hand similar to 
Eq. (|l.6| ), without worrying about its origin. In fact, historically, such an operator was 
introduced suppressed by the Planck scale (rather than the right-handed neutrino mass 
scales) by Weinberg in order to account for small neutrino masses [74]. The problem 
is that such a Planck scale suppressed operator would lead to neutrino masses of order 
10-^eV which are too small to account for the two heavier neutrino masses (though it 
could account for the lightest neutrino mass). To account for the heaviest neutrino mass 
requires a dimension 5 operator suppressed by a mass scale of order 3 x 10^'^ GeV if the 
dimensionless coupling of the operator is of order unity, and the Higgs vacuum expectation 
value (VEV) is equal to that of the Standard Model. 

There are several different kinds of seesaw mechanism in the literature. In this review 
we shall focus on the simplest type I seesaw mechanism, which we shall introduce below. 
However for completeness we shall also discuss the type II and III seesaw mechanisms and 
the double seesaw mechanism, as well as the inverse and linear seesaw mechanisms. 

4.1 Type I seesavi^ 

Before discussing the seesaw mechanism it is worth first reviewing the different types of 
neutrino mass that are possible. So far we have been assuming that neutrino masses are 
Majorana masses of the form 

mlLWi^i , (4.2) 

where vl is a left-handed neutrino field and f£ is the CP conjugate of a left-handed neutrino 
field, in other words a right-handed antineutrino field. Such Majorana masses are possible 
since both the neutrino and the antineutrino are electrically neutral and so Majorana 
masses are not forbidden by electric charge conservation. For this reason a Majorana mass 
for the electron would be strictly forbidden. However such Majorana neutrino masses 
violate lepton number conservation. Assuming only Higgs doublets, they are forbidden 
in the Standard Model at the renormalisable level by gauge invariance. The idea of the 
simplest version of the seesaw mechanism is to assume that such terms are zero to begin 
with, but are generated effectively, after right-handed neutrinos are introduced [33]. 

If we introduce right-handed neutrino fields then there are two sorts of additional 
neutrino mass terms that are possible. There are additional Majorana masses of the form 

MRR^iyR , (4.3) 

where is a right-handed neutrino field and i/^ is the CP conjugate of a right-handed 
neutrino field, in other words a left-handed antineutrino field. In addition there are Dirac 
masses of the form 

miLRWi^R ■ (4.4) 
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Figure 5: Diagram illustrating the type I seesaw niechanism. 



Such Dirac mass terms conserve lepton number, and are not forbidden by electric charge 
conservation even for the charged leptons and quarks. 

Once this is done then the types of neutrino mass discussed in Eqs. ( 0| , 47^ (but not 



Eq. (|4.2| ) since we assume no Higgs triplets) are permitted, and we have the mass matrix 



Since the right-handed neutrinos are electroweak singlets the Majorana masses of the right- 
handed neutrinos MpiR may be orders of magnitude larger than the electroweak scale. In 
the approximation that Mrji ^ ttilr the matrix in Eq. ( [4. 51 ) may be diagonalised to yield 
effective Majorana masses of the type in Eq. (^), 

"^LL = -mLRM^lm^R . (4.6) 

The effective left-handed Majorana masses mJ^^j^ are naturally suppressed by the heavy 
scale MjiFi. In a one family example if we take miR ^ My^r (where M\y is the mass of 
the W boson) and Mrr ~ Mgut then we find w.^^ ~ 10"'^ eV which looks good for 
solar neutrinos. Atmospheric neutrino masses would require a right-handed neutrino with 
a mass below the GUT scale. 

The type I seesaw mechanism is illustrated diagrammatically in Fig. |5[ It can be 
formally derived from the following Lagrangian 

C = -ITZUILRVR - ^^VrMrrvr + H.c. , (4.7) 

where vl represents left-handed neutrino fields (arising from electroweak doublets), vr rep- 
resents right-handed neutrino fields (arising from electroweak singlets), in a matrix notation 
where the m^/j matrix elements are typically of order the charged lepton masses, while the 
MjiR matrix elements may be much larger than the electroweak scale, and maybe up to 
the Planck scale. The number of right-handed neutrinos is not fixed, but the number of 
left-handed neutrinos is equal to three. Below the mass scale of the right-handed neutrinos 
we can integrate them out using the equations of motion 

!^ = 0, (4.8) 

VR 

which gives 

^R = -^mLRMj^}^ , un = -M^^lfn^j^^ . (4.9) 
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Substituting back into the original Lagrangian we find 

C = -^171^11^1^1+ H.c. , (4.10) 

with m'^j^ as in Eq. ([1.6|). 

4.2 Minimal seesaw extension of the Standard Model 

We now briefly discuss what the Standard Model looks like, assuming a minimal seesaw 
extension. In the Standard Model Dirac mass terms for charged leptons and quarks are 
generated from Yukawa type couplings to Higgs doublets whose vacuum expectation value 
gives the Dirac mass term. Neutrino masses are zero in the Standard Model because right- 
handed neutrinos are not present, and also because the Majorana mass terms in Eq. ( [4.2| ) 
require Higgs triplets in order to be generated at the renormalisable level. The simplest 
way to generate neutrino masses from a renormalisable theory is to introduce right-handed 
neutrinos, as in the type I seesaw mechanism, which we assume here. The Lagrangian 
for the lepton sector of the Standard Model containing three right-handed neutrinos with 
heavy Majorana masses is^ 



+ H.C., (4.11) 



where Ca/s = —(-pa^ ei2 = 1; and the remaining notation is standard except that the 3 right- 
handed neutrinos have been replaced by their CP conjugates z/?, and M^^j^ is a complex 
symmetric Majorana matrix. When the two Higgs doublets get their VEVs < H^ >= Vu, 
< H\ >= Vd, where the ratio of VEVs is defined to be tan/3 = Vu/vd, we find the terms 

Cn.^s = -VdX^e^f^ - VuYi^iy.u^ - iM'iRiyii^j + H.c. . (4.12) 
Replacing CP conjugate fields we can write in a matrix notation 

>Cmass = -eLVdY*eR - 17ZVuY*VR - ^J^rMrrI'r + H.c. . (4.13) 

It is convenient to work in the diagonal charged lepton basis 

diag{me, m^, rur) = Ve^VdY*V^^ , (4.14) 

and the diagonal right-handed neutrino basis 

diag(Mi, M2, Ms) = V,^M*j,j,Vj^ , (4.15) 

where 14^,V^^,V^^ are unitary transformations. In this basis the neutrino Yukawa cou- 
plings are given by 

Y, = Ve,Y:vJ^ , (4.16) 



^We introduce two Higgs doublets to pave the way for the supersymmetric Standard Model. For the same 
reason we express the Standard Model Lagrangian in terms of left-handed fields, replacing right-handed 
fields tpR by their CP conjugates tp'^, where the subscript R has been dropped. In the case of the minimal 



standard seesaw model with only one Higgs doublet, the other Higgs doublet in Eq. (4.11) is obtained by 
charge conjugation, i.e. Hd = H^. 
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and the Lagrangian in this basis is 



- {I'eL T^iiL ^tl)VuYu{vri yR2 l^Rs) 

- ^(J^m '^R2 i^R3)diag{Mi, M2, M3){uRi ur2 VR^f + H.c. . (4.17) 
After integrating out the right-handed neutrinos (the seesaw mechanism) we find 

>Cmass =- ieZJiZTZ)diag{me,m^,mr){eRfj.RTR)'^ 

- ^{V^ VJZ ^)rn'iL{i^el J^f^l i^rif + H.c. , (4.18) 

where the hght effective left-handed Majorana neutrino mass matrix in the above basis is 
given by the following seesaw formula which is equivalent to Eq. (4.6), 

<L = -vlY,diag{M^\M^\M^')Yj . (4.19) 

In this basis the type I seesaw mechanism reproduces the dimension 5 operator in Eq. ( [4.1| ) 
with 

K = Yu diag(Mf \ M2"\ Mg-^) yJ . (4.20) 
4.3 Sequential right-handed neutrino dominance 

In this subsection we show how the type I seesaw mechanism may lead to a neutrino mass 
hierarchy with large solar and atmospheric mixing angles, and a reactor angle of order 
1112/ ni'i via a simple and natural mechanism known as sequential dominance (SD).® First 
consider the case of single right-handed neutrino dominance where only one right-handed 
neutrino f| of heavy Majorana mass M3 is present in the seesaw mechanism, namely the 
one responsible for the atmospheric neutrino mass [51,52]. We work in the basis of the 
previous subsection where the right-handed neutrinos and charged lepton mass matrices 
are diagonal. If the single right-handed neutrino couples to the three lepton doublets Li in 
the diagonal charged lepton mass basis as, 

Hu{dLe + eL^ + fLr)vl (4.21) 

where d,e,f are Yukawa couplings (assumed real for simplicity^), and H^ is the Higgs 
doublet, where it is assumed that d <^ e, f. so that the seesaw mechanism yields the 
atmospheric neutrino mass, 

"^3«(e2 + /')^, (4.22) 

where Vu = {Hu). Then the reactor and atmospheric angles are approximately given by 
simple ratios of Yukawa couplings [51,52], 

d „ e 
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tan6'23~-. (4.23) 



®SD should not be confused with an alternative mechanism proposed by Smirnov which is based on the 
premise that there are no large mixing angles in the Yukawa sector, and does not discuss any mechanism 
for achieving this involving right-handed neutrinos [75]. 

'^The full results including phases are discussed in [51,52] and summarised for the case of d = in [53]. 
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According to SD [51, 52] the solar neutrino mass and mixing are accounted for by 
introducing a second right-handed neutrino with mass M2 which couples to the three 
lepton doublets Lj in the diagonal charged lepton mass basis as, 

Hu{aLe + bL^ + cLrM, (4.24) 

where a,b,c are Yukawa couplings (assumed real for simplicity). Then the second right- 
handed neutrino is mainly responsible for the solar neutrino mass, providing 

(a,6,c)VM2«(e,/)VM3, (4.25) 

which is the basic SD condition. Assuming this, then the seesaw mechanism leads to the 
solar neutrino mass, 

ma ~ (a2 + (C236 - S23cf) ^ , (4.26) 

and the solar neutrino mixing is approximately given by a simple ratios of Yukawa cou- 
plings [51,52], 

tan0i2«7 r. (4.27) 

There is an additional contribution to the reactor angle of the form [51,52], 

There may also be a third right-handed neutrino but with completely subdominant contri- 
butions to the seesaw mechanism, and hence it may ignored to leading order. ^ 

Let us summarise what SD achieves. With the assumption in Eq. ( [4.25|) , SD predicts 
a neutrino mass hierarchy, together with solar and atmospheric mixing angles which are 
independent of neutrino mass. Since they only involve ratios of Yukawa couplings they 
may easily be large. On the other hand the reactor angle has two contributions, one 
proportional to a ratio of Yukawa couplings which may be small if d <C e, while the other 
one gives a contribution O (771-2 /?ti-3) which is by itself of the correct magnitude to account 
for the reactor angle (even if d = 0). The origin of these conditions and assumptions may 
be due to family symmetry as we will discuss. 

4.4 Other seesaw mechanisms 

One might also wonder if the seesaw mechanism with right-handed neutrinos is the only 
possibility? In fact it is possible to generate the dimension 5 operator in Eq. (|1.6|) by 
the exchange of heavy Higgs triplets of SU{2)l, referred to as the type II seesaw mech- 
anism [34] or the exchange of heavy SU{2)l triplet fermions, referred to as the type III 
seesaw mechanism [35]. 

In the type II seesaw the general neutrino mass matrix is given by 



*The contributions of the third sub-subdominant right-handed neutrino to the mixing angles has been 
considered in [76]. 
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Figure 6: Diagram leading to a type II contribution to the neutrino mass matrix via an 

induced VEV of the neutral component of a triplet Higgs A. 



Under the assumption that the mass eigenvalues Mj of M^r are very large compared 
to the components of and uilr, the mass matrix can approximately be diagonalised 
yielding effective Majorana masses 

<L ^ rnii + mi^ , (4.30) 

with 

= -rriLR M'}^ mlji , (4.31) 

for the light neutrinos. The direct mass term can also provide a naturally small 

contribution to the light neutrino masses if it stems e.g. from a seesaw suppressed induced 
VEV, see Fig. ^. The general case, where both possibilities are allowed, is also referred to 
as the type II seesaw mechanism. 

Alternatively the seesaw can be implemented in a two-stage process by introducing 
additional neutrino singlets beyond the three right-handed neutrinos that we have consid- 
ered so far. It is useful to distinguish between "right-handed neutrinos" which carry 
B — L and perhaps form SU{2)r doublets with right-handed charged leptons, and "neutrino 
singlets" S which have no Yukawa couplings to the left-handed neutrinos, but which may 
couple to vpi. If the singlets have Majorana masses Mss, but the right-handed neutrinos 
have a zero Majorana mass Mr^ = 0, the seesaw mechanism may proceed via mass cou- 
plings of singlets to right-handed neutrinos Mas- In the basis {vl, vr, S) the mass matrix 
is 

( ruLR \ 
ml^ Mrs ■ (4.32) 

There are two different cases often considered: 

(i) Assuming Mss ^ Mrs the light physical left-handed Majorana neutrino masses 
are then, 

"T-LL = rriLRM-j^mlji , (4.33) 

where 

Mrr = MrsM^ImIs. (4.34) 
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This is called the double seesaw mechanism [40] . It is often used in GUT or string inspired 
neutrino mass models to explain why Mrji is below the GUT or string scale. 

(ii) Assuming Mss «C Mrs, the matrix has one light and two heavy quasi-degenerate 
states for each generation. The mass matrix of the light physical left-handed Majorana 
neutrino masses is, 



which has a double suppression. 

In the limit that Mss ~^ all neutrinos become massless and lepton number sym- 
metry is restored. Close to this limit one may have acceptable light neutrino masses for 
Mrs ~ 1 TeV, allowing a testable low energy seesaw mechanism referred to as the inverse 



one obtains another version of the low energy seesaw mechanism called the linear seesaw 
mechanism. 



5. Finite group theory in a nutshell 

In this section, we give a first introduction to finite group theory, using the permutation 
group of three objects 5*3 as an example, and later generalising the discussion to include 
all finite groups with triplet representations. Readers who are familiar with finite group 
theory may wish to skip this section. 

5.1 Group multiplication table 

Non-Abelian discrete symmetries appear to play an important role in understanding the 
physics of flavour. In order for this pedagogical review to be self-contained, we give a 
brief introduction into the main mathematical concepts of finite group theory. Many more 
details can be found, for instance, in the recent textbook by Ramond [78] which provides 
an excellent survey of the topic particularly aimed at physicists. 

Finite groups G consist of a finite number of elements g together with a binary op- 
eration that maps two elements onto one element of G. In the following we use the term 
multiplication for such an operation. By definition, a group must include an identity ele- 
ment e as well as the inverse of a given element g. Furthermore, the multiplication must 
be associative, meaning that the product of three elements satisfies (91(72)53 = 51(9253)- 
Groups are called Abelian if gig2 = 5251 for all group elements, while the elements of 
non-Abelian groups do not satisfy this trivial commutation relation in general. We shall 
only be interested in non-Abelian groups from now on. 

The most basic way of defining a group is given in terms of the multiplication table, 
where the result of each product of two elements is listed. In the case of the smallest 



<L = mLRM]^s^MssM;^^mlR, 



(4.35) 



seesaw mechanism. If one allows the 1-3 elements 
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non-Abelian finite group, the permutation group 53, we have: 
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The six elements are classified into the identity clement e, elements bi whose square is e 
and finally elements Oj for which the square docs not yield e but, as can be seen easily, 
the cube does. It is generally true for any finite group that there exists some exponent n 
for each element g such that g'" = e. The smallest exponent for which this holds is called 
the order of the element g. This is not to be confused with the order of a group G which 
simply means the number of elements contained in G. 

5.2 Group presentation 

Clearly, the definition of a finite group in terms of its multiplication table becomes cum- 
bersome very quickly with increasing order of G. It is therefore necessary to find a more 
compact way of defining G. Noticing that all six elements of 5*3 can be obtained by mul- 
tiplying only a subset of all elements, we arrive at the notion of generators of a group. 
Denoting ai = a and bi = b, we obtain 02 = as well as 62 = ab and 63 = ba. In other 
words, a and b generate the group S3. Being the group of permutations on three objects 
which is isomorphic to the group of symmetry transformations of an equilateral triangle, a 
corresponds to a 120° rotation and 6 to a reflection. This observation leads to the deflnition 
of S3 using the so-called presentation 

< a, b \ = b'^ = e, bab~^ = > , (5.1) 

where the generators have to respect the rules listed on the right. Depending the these 
presentation rules, a group can be defined uniquely in a compact way. Unfortunately, 
such an abstract definition of a group is not very useful for physical applications as it 
does not show the possible irreducible representations of the group. We therefore quickly 
continue our journey through the fields of finite group theory towards the important notion 
of character tables. 

5.3 Character table 

In order to understand the meaning of a character table, is it mandatory to introduce the 

idea of conjugacy classes and irreducible representations. Conjugacy classes are subsets 
of elements of G which are obtained from collecting all those elements related to a given 
element gi by conjugation ggig~^, for all g & G. The union of all possible conjugacy classes 
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is nothing but the set of ah elements of G. In the case of 6*3 we find three different classes, 



1C\1) = {glg-'lgeSs} = {1} , 
2C3(a) = {gag-'\geS3} = {a, a^} , 

3C\b) = {gbg-' U G 53} = {b, ab, ba} . (5.2) 



Here we have used the notation NiC^''{gi), where gi is an element of the class, Ni gives 
the number of different elements contained in that class, and Ui denotes the order of these 
elements, which is identical for all ggig^^ with g £ G. 

The other ingredient for constructing a character table is the set of possible irreducible 
representations of the group G. In general non-Abelian groups can be realised in terms of 
r X r matrices, where the positive integers r depend on the group. Then, the abstract gener- 
ators of a group are promoted to matrices which satisfy the presentation rules. Such matrix 
representations are called reducible if there exists a basis in which the r x r matrices of all 
generators of G can be brought into the same block diagonal form. If this is not possible, 
the representation is called irreducible. Clearly, the trivial singlet representation 1, where 
all generators of G are identically 1, satisfies any presentation rule and is thus an irreducible 
representation of all groups. This trivial example shows that irreducible representations do 
not necessarily have to be faithful, i.e. multiplying the matrices corresponding to the group 
generators can give a smaller number of different matrices than the order of G. In the case 
of 5*3, the irreducible representations compatible with the presentation rules of Eq. ( ^Tl] ) 
take the form 



1 : a = 1 , 
1' : a = 1 , 



-1 , 



2 : a 



/ 2TTi 

e— 

2Tvi 

, e~ — 




(5.3) 



The fact that 5*3 has three irreducible representations and also three conjugacy classes is 
not a coincidence. It is generally true that the number of irreducible representations of a 
finite group is equal to the number of its conjugacy classes. Moreover, summing up the 
squares of the dimensions of all irreducible representations always yields the order of the 
group G. For example, in S3 we get 1^ + 1^ + 2^ = 6. These two facts can be used to work 
out all irreducible representations of a given group G. 

In the case of irreducible representations r with r > 1, the explicit matrix form of 
the generators depends on the basis. In order to obtain a basis independent quantity, one 
defines the character Xgt of the matrix representation of a group element g^ to be its trace. 
Since the elements of a conjugacy class are all related by ggig^^ with (7 € G, it is meaningful 
to speak of the character Xi of the elements of a conjugacy class i. Therefore one can 
define the (quadratic) character table where the rows list the irreducible representations 
and the columns show the conjugacy classes. Using Eq. ( |5.3| ), we easily find the following 
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character table of S3. 



S3 


1C\1) 


2C^a) 


3C2(6) 




1 


1 


1 


xP'' 


1 


1 


-1 


xf 


2 


-1 






Defining a group in terms of its character table is much more suitable for physical applica- 
tions than the previous two definitions. First, it immediately lists all possible irreducible 
representations which might be used in constructing particle physics models. Secondly, 
it is also straightforward to extract the Kronecker products of a finite group G from its 
character table. 

5.4 Kronecker products and Clebsch-Gordan coefficients 

Multiplying arbitrary irreducible representations r and s 

r(g)s = ^ti(r,s,t)t, (5.4) 
t 

one can calculate the multiplicity d{r, s, t) with which the irreducible representation t 
occurs in the product by 

d(r,s,t) = lj2^i-xM\f , (5.5) 

i 

where the sum is over all classes. N denotes the order of the group G and the asterisk 
indicates complex conjugation. With this, we obtain the following non-trivial Kronecker 
products from the S3 character table, 

I'^l' = 1 , 
1' ® 2 = 2 , 
2(g)2 =1 + 1' + 2 . 

The Kronecker products are necessarily independent of the bases of the irreducible rep- 
resentations r with r > 1. When formulating and spelling out the details of a model, 
particular bases have to be chosen by hand. With the bases fixed, it is possible to work out 
the basis dependent Clebsch-Gordan coefficients of a group. Denoting the components of 
the two multiplet of a product by cij and Pj , the resulting representation with components 
7fe are obtained from, 

lk = Y^ Cij ai I3j , (5.6) 

where arc the Clebsch-Gordan coefficients. These are determined by the required trans- 
formation properties of the components 7^ under the group generators. In the case of S3, 
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using the basis of Eq. (|5.3D , one gets, 

I'Ol' 1 

l'(g)2 2 

2 2^1 
2 02 1' 

2 2 2 

where ai refers to the first factor of the Kronecker product and f3j to the second. We 
conclude our discussion of the most important concepts in finite group theory by pointing 
out that - due to the choice of convenient bases - a representation which is real (that is 
for which there exists a basis where all generators are explicitly real) may have complex 
generators. This is for instance the case for the doublet of ^3 in the basis of Eq. ( |5.3| ). 

5.5 Finite groups with triplet representations 

For applications in flavour physics, we are interested in finite groups with triplet representa- 
tions. They can be found among the subgroups of SU{3) and fall into four classes [44,79]:^ 

• Groups of the type (Z„ x Zm) x 5*3 

• Groups of the type {Zn x Zm) x Z3 

• The simple groups and PSL2{7) [81] plus a few more "exceptional" groups [82] 

• The double covers of the tetrahedral (^4), octahedral (84,) and icosahedral (^5) 
groups 

The latter are subgroups of SU{2), whose triplet representations are identical to the triplets 
of the respective rotation groups (which in turn are already included in the other classes). 
Many of the physically useful symmetries are special cases within these general classes. 
For instance, ^4, the natural symmetry of tri-bimaximal mixing in direct models, see 
Subsection is isomorphic to A{6n?) = (Z„ x Z„) x 5*3 with n = 2. The presentation 
rules of A(6n^) can be given in terms of four^*^ generators, a,b,c,d [83], 

a3 = 52 = (a6)2 = c" = d'^ = 1 , cd = dc , 

aca'^ = c^^d^^ , ada~^ = c , bcb^^ = d^^ , bdb^^ = c^^ . (5-7) 

^Subgroups of (7(3) can be derived from SU{3) [44], however, a complete classification is still lacking [80]. 
^"in principle, the presentation can be easily formulated with only three generators by expressing either 
c or d in terms of the other three generators. 
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ai/32 + a2/5i , 
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The dimensions of all irreducible representations can only take values 1, 2, 3 or 6. A faithful 
triplet representation is found, e.g., in the following set of matrices [83], 









/o 1^ 




(v 


o\ 




/l \ 


a = 


1 




10 




T]-^ 





, d = 


7? 




[lOOj 




[lOOj 




i^O 






^0 7]-^) 



(5.8) 



where rj = e~ . With n = 2 this triplet representation is explicitly real, and therefore 
does not correspond to the basis in which the 5*4 order three generator T is diagonal and 
complex, cf. Section ^. To make connection to the 54 triplet generators S, U and T as 
listen in Appendix |^, we have to perform the basis transformation [84] , 



S = wdw 



U = w (aba ^) w 



T 



w aw 



where 



w 



1 

7i 



/i 



1 1 



1 w w 
1 w 



with 



2-Ki 

e 3 



(5.9) 



(5.10) 



This shows how the tri-bimaximal Klein symmetry Z2 x Z2 of the neutrino mass matrix 
in the diagonal charged lepton basis, generated by S and U of Eq. ( |6.8D , is inherited from 
A(24) = {Z2 X Z2) X 53: one Z2 factor (namely S) originates from the first factor, Z2 x Z2, 
and the other (namely U) is derived from the second, 53. We remark in passing that the 
smallest group within the series A(6n^) containing sextets is A(96). 

Another series of groups can be obtained from the presentation of Eq. ( |5.7D by simply 
dropping the generator 6, and consequently all conditions involving h [85]. This results 
in the groups A(3n^) = (Z„ x Z„) x Z3 which only allow for irreducible representations 
of dimension 1 and 3. The case with n = 2 generates the tetrahedral group A4, and the 
faithful triplet representation is the same as in the case of only without the h ox U 

3 we obtain the group A (27) which has also been 



generator, cf. Eqs. (5^,0|). With n 
applied successfully as a family symmetry in indirect models [86-89]. 

A third series is obtained from the second class of groups, (Z„ x Zm) x Z3, by setting 
m = 1. The presentation of this series of groups T„ = Z„ x Z3 reads [90] 



1 



aca 



(5.11) 



where the integer k must satisfy 1 + A; + /c^ = mod n. One can easily check that, with 
r] = e~ , a faithful triplet representation is given by 



(5.12) 



Popular examples of such groups include Ty [91] and T13 [92], both of which do not include 
a Z2 subgroup so that the Klein symmetry of the neutrino mass matrix cannot be obtained 
from these groups in a direct or semi-direct way, see Subsection Yet, from the model 
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Figure 7: Examples of subgroups of SU{3) with triplet representations discussed in this review. 
A line connecting two groups indicates that the smaller is a subgroup of the bigger one. 

building point of view it can still be useful to change to a basis in which the order three 
generator becomes diagonal [93], analogously to the case of S4. In Appendix |p we list the 
generators and Clebsch-Gordan coefficients of the groups S4, A4 and T7 in the T diago- 
nal basis. Their relation to SU{3) and some of its subgroups is schematically illustrated 
in Fig. 0. 



6. Discrete family symmetries and model building approaches 
6.1 Family symmetries and flavons 

The masses and mixings of the three families of quarks and leptons result from the form 
of the respective Yukawa matrices formulated in the flavour basis. Is there an organising 
principle which dictates the family structure of these Yukawa couplings? While this review 
takes the view that the observed mass and mixing patterns can be traced back to a family 
symmetry, we remark that some authors answer this question negatively, referring to a 
landscape of parameter choices out of which Nature has picked one that is compatible with 
the experimental measurements. In particular, the observation of a large reactor angle 
has been interpreted as a sign for an anarchical neutrino mass matrix [94]. Following 
the symmetry approach, it is clear that the family symmetry must be broken in order 
to generate the observed non-trivial structures. This is achieved by means of Higgs-type 
fields. These so-called flavon fields <j) are neutral under the SM gauge group and break 
the family symmetry spontaneously by acquiring a VEV. This VEV in turn introduces an 
expansion parameter 

A ' 



(6.1) 
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(A denotes a high energy mass scale) which can be used to derive hierarchical Yukawa 
matrices, possibly with texture zeroes. 

Family symmetries are sometimes also called horizontal symmetries, as opposed to 
GUT symmetries which unify different members within a family. It is possible to impose 
Abelian or non-Abelian family symmetries. The former choice goes back to the old idea 
of Froggatt and Nielsen [95] to explain the hierarchies of the quark masses and mixings in 
terms of an underlying {7(1)fn symmetry. In such a framework, the three generations of 
(left- and right-handed) quark fields carry different charges under C/(1)fn such that the 
usual Yukawa terms have positive integer charges. Depending on the involved generations, 
this can be compensated by multiplying nij factors of the flavon field (f) which conventionally 
has a C/(1)fn charge assignment of —1. The resulting terms which give rise to the usual 
Yukawa interactions then take the form 



^ I qLiqRi H , (6.2) 

where i and j are family indices and Cij denote undetermined order one coefficients. In- 
serting the flavon VEV then generates the Yukawa couplings Yij = Cije"''^ which become 
hierarchical if the C/(1)fn charges are chosen appropriately. We emphasise that this ap- 
proach is only useful for explaining hierarchical structures as the order one coefficients Cij 
remain unspecified. 

In order to accurately describe non-hierarchical family structures such as the observed 
peculiar lepton mixing pattern, it is necessary to impose a non-Abelian family symmetry. 
The three generations of quarks and leptons can then be unified into suitable multiplets 
(i.e. irreducible representations) of the family symmetry G. An intimate connection of 
all three families is provided if G contains triplet representations, ^ = (^i, V'2; ^3)"^ ~ 3. 
Requiring irreducible triplet representations, the possible choices for G are limited to U (3) 
and subgroups thereof. To illustrate the idea, we sketch the essential steps using the 
example of SU{3) [96], applied to the Weinberg operator HuL^ kLHu, cf. Eq. (|l.6|) . The 
three generations of lepton doublets are unified into a triplet of SU (3) while the Higgs Hu 
is assumed to be a singlet. In order to construct an SU{3) invariant operator, a flavon field 
(j) transforming in the 3 of SU (3) can be introduced, leading to the term 

HuL^(t)(f^LHu . (6.3) 

The VEV of the flavon field (j) will now correspond to a vector with a particular alignment, 
i.e. {(j)) oc {a,b,c)'^, where a,b,c take numerical values dictated by the scalar potential. 
Inserting this vacuum configuration into the factor (fxj)^ of Eq. ( |6.3| ) will generate a contri- 
bution to the neutrino mass matrix which is proportional to 

^a^ ah ac^ 

ha 6^ he . (6.4) 
\ca ch (? j 

Assuming simple fiavon alignments, it is possible to relate all entries of the neutrino mass 
matrix in a particular way so that special mixing patterns can be explained. In practice, at 
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least two flavons with different ahgnments have to be imposed in order to avoid degeneracies 
in the neutrino masses. Although it is possible to obtain simple and predictive flavon 
alignments in models based on continuous family symmetries such as SU{3) and 50(3), 
the problem of vacuum alignment can be solved in a significantly simpler and more natural 
way by imposing a discrete family symmetry instead. In the following we therefore focus 
our attention on non-Abelian discrete family symmetries with triplets. 

6.2 The Klein symmetry of the neutrino mass matrix 

The PMNS mixing is dictated by the structure of charged and neutral lepton mass ma- 
trices in a weak eigenstate basis. More precisely it is obtained as the mismatch of the 
transformations on the two left-handed lepton states needed to bring the charged lepton 
and the neutrino mass matrices into diagonal form. In order to easily reach a physical 
interpretation, it is convenient to work in a basis in which the charged leptons are diag- 
onal, or approximately diagonal. The latter is useful in GUT model building where the 
non-diagonal hierarchical down-type quark mass matrix, required for the observed CKM 
mixing, is directly related to the charged lepton mass matrix which, as a consequence, is also 
not completely diagonal. The total PMNS mixing will then be predominantly determined 
by the neutrino mass matrix, and small charged lepton corrections, see Subsection |3.4| , 
will have to be taken into account separately, leading to characteristic mixing sum rules as 
explained in Subsection |3.5| . 

With this in mind, one can hope to obtain clues on the nature of the underlying 
family symmetry by studying the symmetry properties of the neutrino mass matrix in the 
basis of (approximately) diagonal charged leptons. Assuming neutrinos to be Majorana 
rather than Dirac particles, their mass matrix is always symmetric under a Klein symmetry 
Z2 X Z2- This follows from the obvious observation that the diagonalised neutrino mass 
matrix m'^'^'^^ is left invariant under the transformation 

/(-1)P \ 

(-!)<? , (6.5) 

\ {-i)p+y 

where p and q take the integer values and 1. The explicit form of the Klein symmetry of 
the non-diagonalised neutrino mass matrix m'^j^, expressed in terms of 3 x 3 matrices, can 
then be determined as 

Kp^q = C^PMNS-^P,g^PMNS ' (6-6) 

where C/pMNS is (approximately) the unitary PMNS mixing matrix. The matrices Kp^q 
form a group of four elements whose squares yield the identity element Kq q. The fact 
that the neutrino mass matrix is symmetric under a transformation by i^p,q can be easily 
verified using Eq. ( 3.14| ), 



^p,q "^LL ^p,q = "^LL ' With Kp^q 



^p,q''^LL ^p,q — ^PMNS^p,qiUpMNS ^LL ^PMNs)^p,q^PMNS 

= ^PMNs(-^p^g "^LL ^ ^p,q)UpMm 

TT i^.diag rrT 

— (^PMNS ^LL '^PMNS 



(6.7) 
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We point out that the Z2 x Z2 symmetry of Eq. ( p.q) exists for any choice of PMNS mixing. 
In the remainder of this review we wih denote the two generators of this Klein symmetry by 
S and \J . Particularly simple forms of these generators are obtained when C/pMNS features 
a simple mixing pattern. For instance, in the case of tri-bimaximal mixing, C/pMNS = f^TB) 
we find 

. /-I 2 2 \ /l 0\ 



S 



-1 2 



-V 



1 
Voioy 



(6.8) 



Such a symmetry of the neutrino mass matrix is only meaningful if the charged leptons 
are (approximately) diagonal. Therefore, it is useful to consider also the (approximate) 
symmetry of the charged lepton mass matrix Mg. As charged leptons are Dirac particles, 
one has to consider the square M^-MX which - if diagonal - is symmetric under a general 
phase transformation T, 



(MsM]) T = MeMl , with T 



/l 



e" 



(6.9) 



\0 
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and m being an integer. The smallest value of m which also enforces a diagonal MgMe 
is m = 3. Choosing m = 2 would leave room for non- vanishing off-diagonal entries in the 
2-3 block of the squared charged lepton mass matrix, which however could be removed by 
imposing a second T-type generator with permuted diagonal entries. 

We conclude the discussion of the symmetries of the mass matrices by emphasising 
that the T symmetry of the charged leptons can hold exactly in models which are only 
concerned about the lepton sector. In GUT models, which additionally include the quarks, 
such a T symmetry is usually only valid approximately. 



6.3 The direct model building approach 

Family symmetry models can be classified according to the origin of the symmetry of the 
neutrino mass matrix. The neutrino Klein symmetry can arise as a residual symmetry of 
the underlying family symmetry G, in other words, the four elements Kp^g of Eq. ( |6.(^ ) 
are also elements of the imposed family symmetry. Models of this type are called direct 
models [97]. 

In such models, the neutrino mass term involves flavon fields (p'^ whose vacuum align- 
ments break the family symmetry G down to the remnant Klein symmetry of Eq. (|6.6|). 
Schematically this can be expressed as 



A2 



LHuLHu 



with S{(t)l = U{(l)'') 



(6.10) 



where the flavon enters only linearly. Depending on the family symmetry, there are several 
neutrino-type flavons 4>'^ which contribute to the neutrino mass matrix and furnish different 
representations of G. Therefore, S and U stand for the respective Klein generators in the 
representation of (j)'^. In particular, if the group theory allows LHuLHu to be contracted 
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Figure 8: A sketch of the direct model building approach. The charged lepton sector is (approxi- 
mately) diagonal either due to a remnant (approximate) T symmetry or simply by construction. 

to a singlet 1 of G, then it is possible to introduce a flavon in the 1 of G with S = U = 1. 
Such a trivial flavon would clearly not break the Klein symmetry, in fact, it does not even 
break G. The same Klein symmetry can be realised straightforwardly in direct models 
based on the type I seesaw mechanism. There, the Dirac neutrino term does not involve 
a flavon field, and therefore does not break G at all, whereas the Majorana neutrino mass 
term involves the S and U preserving flavons linearly, 

C ^ Lv^Hu + <t)''v^y'' . (6.11) 

Application of the type I seesaw formula yields an effective light neutrino mass matrix 
which is again symmetric under S and U . Analogously, the charged lepton sector often 
involves a flavon (j)^ which breaks G (approximately) to the symmetry generated by T, 

~ ^LtHd , with r(/) « (/) . (6.12) 

The direct approach is schematically illustrated in Fig. However, we remark that accord- 
ing to our classification, the charged lepton sector may be diagonal simply by construction. 

Assuming the T symmetry to be exact, the class of direct models clearly requires G 
to contain both the Klein symmetry of the neutrino sector as well as the symmetry of the 
charged leptons. The minimal symmetry group for which this is satisfied can be determined 
by simply calculating all possible products of the matrices S, U and T. In the tri-bimaximal 
case, i.e. with the generators of Eq. ( |6.8| ) for the neutrino Klein symmetry, and a charged 
lepton symmetry T of Eq. ( |6.9D with m = 3 one obtains a total of 24 different matrices 
which form a finite group isomorphic to 5*4, the group of permutation on four objects. It is 
interesting to note that a different choice of m does not necessarily yield a finite group. In 
fact, with S and U of Eq. (^) one can easily show, using the computer algebra programme 
GAP [98], that no finite group is obtained for reasonably small m ^ 3. We have checked 
the validity of this statement for m < 30, which is even true for the case with m = 2. In 
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that sense, 5*4 is the natural symmetry group of tri-bimaximal mixing in direct models. Of 
course, any group that contains 5*4 as a subgroup could be applied equally well. 

In addition to the pure class of direct models, there are semi-direct models in which one 
Z2 factor of the Klein symmetry can be identified as a subgroup of G, while the other Z2 
factor arises accidentally. The Savons of semi-direct models appear linearly in the neutrino 
mass term, similar to Eq. ( |6.10D , and break G down to one of its Z2 subgroups. An example 
of such a model is provided by the famous Altarelli-Feruglio ^4 model of tri-bimaximal 
mixing [30,99] . ^4 is the group of even permutations on four object, and as such a subgroup 
of 54. It can be obtained from ^4 by simply dropping the \J generator. Not being a part of 
the underlying family symmetry, it is therefore evident that the \J symmetry of Eq. ( |6.^ ) 
must arise accidentally. 



6.4 The indirect model building approach 

In the class of indirect models, no Z2 factor of the Klein symmetry of Eq. ( |6.6| ) forms a 
subgroup of G. Models of this class are typically based on the type I seesaw mechanism 
together with the assumption of sequential dominance, see Subsection |4.3| . Here, the main 
role of the family symmetry consists in relating the Yukawa couplings d, e, / of Eq. (4.21) as 
well as a, 6, c of Eq. ( 4.24D by introducing triplet flavon fields which acquire special vacuum 
configurations. The directions of the flavon alignments are determined by the G symmetric 
operators of the flavon potential [97]. 

Working in a basis where both the charged leptons as well as the right-handed neutri- 
nos are diagonal, the leptonic flavour structure is encoded in the Dirac neutrino Yukawa 
operator. The triplet flavons of indirect models enter linearly in this term, 



i 



(6.13) 



where A is a cut-off scale and the sum is over the number of right-handed neutrinos. Adopt- 



ing the notation of Subsection O, extended to include a third right-handed neutrino z^J, 
we obtain the Dirac neutrino Yukawa matrix by inserting the flavon VEVs into Eq. ( 6.13| ) . 
Suppressing the dimensionless couplings of the Dirac neutrino terms for notational clarity, 
we get 
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(6.14) 



The columns of the Dirac neutrino Yukawa matrix are therefore proportional to the vacuum 
alignments of the flavons fields (p^. The effective Majorana operators of the light neutrinos 
can be derived from this using the seesaw formula of Eq. (|4.6| ), yielding 
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eff 



^ m j_ 
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(6.15) 



Note that the flavons enter the effective neutrino mass terms quadratically. In models with 
sequential dominance, the three contributions to the effective light neutrino mass matrix 
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Figure 9: A sketch of the indirect model building approach. The charged lepton sector is (approx- 
imately) diagonal by construction. 

are hierarchical, and it is often possible to ignore one term (e.g. i = 1) such that the sum 
contains only one dominant and one subdominant contribution. In the class of indirect 
models, the PMNS mixing pattern thus becomes a question of the alignment vectors {<p\)- 
For instance, a neutrino mass matrix that gives rise to tri-bimaximal mixing can be obtained 
using the flavon alignments 
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A V2 



/o 



V- 



.16) 



Note that the resulting columns of the Dirac neutrino Yukawa matrix are proportional 
to the columns of the unitary (in the present case tri-bimaximal) mixing matrix. Such a 
property of the Dirac neutrino Yukawa matrix is generally called form dominance [100].^^ 
Furthermore these alignments are left invariant under the action of the S and U generators 
of Eq. ( |6.8D , up to an irrelevant sign which drops out due to the quadratic appearance 
of each flavon in Eq. ( |6.15| ). Since the family symmetry G does not contain the neutrino 
Klein symmetry, its primary role is then to explain the origin of these or similarly simple 
flavon alignments. We schematically illustrated the indirect approach in Fig. |^. 

For completeness we also mention that there are a few special models not based on 
the type I seesaw mechanism, in which the Klein symmetry of the neutrino sector is not a 
subgroup of the family symmetry, see e.g. [88,93]. In such models, the flavon enters linearly 
in the neutrino mass term, and the Klein symmetry arises accidentally from a combination 
of the Clebsch-Gordan coefficients and suitable flavon alignments of (j)'^. 

6.5 Comments on the classification into direct and indirect models 

Before continuing to discuss the details of the family symmetry breaking, we compare the 
direct approach with the indirect one, see Figs. and U. This classification is purely based 



Exact form dominance entails vanishing leptogenesis [101]. 
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on the origin of the Z2 x Z2 Klein symmetry of the neutrino sector, formulated in a basis 
of (approximately) diagonal charged leptons. In direct models, this symmetry, generated 
by the order two elements S and [/, arises as a subgroup of G, whereas this is not the 
case for indirect models. In both approaches, the family symmetry G has to be broken 
spontaneously by flavon fields acquiring a VEV. The flavon vacuum configuration of direct 
models is dictated by the requirement that S and \J be preserved. In indirect models 
based on the type I seesaw, the vacuum alignment of the flavons enters in the columns of 
the Dirac neutrino Yukawa matrix, thereby generating contributions to the effective light 
neutrino mass matrix of the form proportional to . 

We emphasise that the charged leptons have to be (approximately) diagonal by con- 
struction for the purpose of this classification. In the framework of direct models, this can 
be enforced by demanding a subgroup of G, generated by T, to be (approximately) pre- 
served in the charged lepton sector. However, in grand unified models such a T symmetry 
cannot be exact as it would then apply to the quarks as well. This in turn would entail a 
phenomenologically unacceptable quark sector without CKM mixing. 

7. Direct model building 

In models based on a family symmetry G, the Dirac Yukawa and the Majorana couplings 
are typically generated dynamically from G invariant operators involving one or more flavon 
fields. In general, these flavons can transform in any of the irreducible representations of G. 
For non-Abelian discrete symmetries, the choice is limited to a flnite set of representations. 
With flavons transforming as multiplets of the family symmetry G, the breaking of G and 
with it the family structure of the Dirac Yukawa and the Majorana couplings crucially de- 
pends on the alignment of the flavon VEVs. In this section we discuss general strategies for 
identifying useful flavons alignments in direct models where the family symmetry G is bro- 
ken to a particular subgroup in the neutrino sector. Furthermore, we give explicit examples 
which illustrate ways of deriving vacuum alignments from flavon potentials. We remark 
that throughout this section we assume a diagonal or approximately diagonal charged lep- 
ton mass matrix which may or may not arise as a result of an (approximately) unbroken 
subgroup of G. 

7.1 Flavon alignments in direct models 

In direct models, flavons enter linearly in the terms of the neutrino Lagrangian as shown 
in Eq. ( |6.1C1| ). These operators are invariant under the full family symmetry G. When G 
is broken spontaneously by the flavon fields acquiring a VEV, the smaller Z2 x Z2 Klein 
symmetry, generated by the order two elements 5, C/ € G, still remains intact. The criteria 
on the vacuum alignment of the involved flavon VEVs can therefore be formulated by the 
condition 

S{ct>n = = ir) . (7.1) 

With flavons generally furnishing different representations r of the family group G, the 
explicit matrix form of S and U clearly differs for different r. For a given representation, it 
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is then straightforward to calculate the form of the alignments which satisfy Eq. . This 
procedure can be repeated for all other representations r, but not all of them will necessarily 



yield a solution to Eq. (7.1), meaning that flavons transforming in such representations 
cannot be adopted in the considered case. 

We discuss this strategy of identifying the structure of the flavon alignments in an ex- 
plicit example for the purpose of illustration. Consider the case of an S4 family symmetry. 
The generators S and U of the tri-bimaximal Klein symmetry are listed in Appendix ^ 
for all five irreducible representations, cf. also Eq. ( |6.8[ ). The S generator of the 1, 1' 
and 2 are all trivial, i.e. the identity element. Therefore any vacuum configuration of 
flavons transforming in these representations will leave invariant the Z2 symmetry associ- 
ated with S. The second Z2 symmetry of the Klein symmetry, generated by U is always 
broken by the VEV of a flavon transforming in the 1' since U = — 1 in this case, while it is 
left intact by a flavon in the 1 of 84^. For the two-dimensional representation, one quickly 
finds that the flavon alignment has to be proportional to (1,1)^ in order not to break the 
U symmetry. Turning to the 3 of 5*4, invariance under U entails a flavon alignment of the 
form (0,1,-1)^. Applying the S transformation on such an alignment yields (0,-1,1)^, 
hence, this alignment does not satisfy Eq. (|7.1| ) as it is not an eigenvector of S with eigen- 
value +1. Finally, we discuss flavons transforming in the 3' representation of S4. The 
most general alignment which is left invariant under the U transformation reads {a,b,b)'^. 
Demanding S{a, b, b)'^ = [a, b, 6)^ fixes a = b, showing that an alignment proportional to 
(1, 1, 1)^ leaves invariant both S and U. Collecting the results of this discussion, we have 
shown that flavons transforming as a 1' and 3 of S4 cannot be used to break the family 
symmetry S4 down to the tri-bimaximal Klein symmetry. On the other hand, flavon flelds 
in the 1, 2 and 3' representations can be adopted, where the latter two have to be aligned 
as 




9^3' 



(7.2) 



in order to leave invariant the Z2 symmetries associated with S and U. Here if denotes the 
overall VEV of a flavon (p. Inserting all three flavons into Eq. (6.10), assuming the lepton 
doublets L to transform in the 3 representation, we end up with a neutrino mass matrix 
which comprises three terms. 
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(7.3) 



Using the matrices S and U of Eq. (|6.8p, one can easily check explicitly that S^m'^^S 



U'^m'^^U = m'^^ as required. Clearly, the alignments of Eq. ( |7.2D depend on the chosen 
basis. In particular the basis of the doublet representation could have been chosen differ- 
ently without affecting the basis of the triplet representation (which we fixed by demanding 
a diagonal T generator). This, however, would also change the Clebsch-Gordan coefficients 
such that the form of the neutrino mass matrix in Eq. ( [7.3[ ) remains unchanged. We em- 
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phasise that the same procedure of identifying the flavon ahgnments of direct models can 
be apphed to arbitrary choices of the Klein symmetry. 



7.2 Vacuum alignment mechanism in direct models 

Having determined the alignments required in a given direct models, the next step is to 
derive them from minimising a flavon potential. In the context of direct models, the most 
popular and perhaps natural approach to tackle the problem of the flavon alignment is 
provided by the so-called F-term alignment mechanism [30,99]. The idea is to couple the 
flavons to so-called driving fields in a super symmetric setup. Like flavons, driving fields are 
neutral under the SM gauge group and transform in general in a non-trivial way under the 
family symmetry G. Introducing a C/(1)r symmetry under which the chiral supermultiplets 
containing the SM fermions carry charge +1, allows to distinguish flavons from driving 
fields by assigning a charge of +2 to the latter while keeping the former neutral. With this 
U{l)ii charge assignment, the driving fields can only appear linearly in the superpotential 
and cannot couple to the SM fermions. The set of superpotential operators involving the 
driving fields Xi is usually referred to as the driving or simply flavon potential VFflavon- 
Assuming that supersymmetry remains unbroken at the scale where the flavons develop 
their VEVs, we can obtain the flavon alignments from the terms of Wflavon by setting the 
F-terms of the driving flelds to zero, i.e. 

-n, = = , (7.4) 

by which the scalar potential is minimised. 

To illustrate the F-term alignment mechanism we give two simple examples based on 
the family symmetry group 6*4. First, consider a driving field X\ and a flavon fleld (j)2 
transforming in the 1 and 2 representations of 54, respectively. Expanding the resulting 
term of the driving superpotential in terms of the component fields cj)2,i we obtain 

Xl(j)24>2 = -'^l (02,102,2 + </'2,2</'2,l) = (/)2,l(/'2,2 • (7.5) 

The F-term condition of Eq. (7.4) then gives rise to the following two solutions. 





(7.6) 

Notice that these two alignments are related by the 5*4 symmetry transformation U, while a 
transformation induced by T does not change the alignment but only the phase of the VEV. 
It is a general feature of any G symmetric theory that one particular solution for the flavon 
alignments will automatically imply a whole set of solutions which are related by symmetry 
transformations. However, the reverse need not be true, i.e. there may be cases in which 
two or more solutions exist which are not related through symmetry transformations. 

As a second example let us consider the alignments of Eq. ( [7.21) . One possible way 
to derive these using the i^'-term alignment mechanism consists in introducing two driving 
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fields, one transforming in the 3 of 5*4, the other in tlie 3' [102]. The corresponding terms 
of the flavon superpotential then read 



where gi are dimensionless coupling constants. Denoting the VEVs of ^g,, 
Cj, bi and a, respectively, the F-term conditions take the form 



6^ and 



(7.7) 

PI by 

(7.8) 
(7.9) 

Restricting to solutions in which all of the three flavons develop a VEV, Eq. ( |7.8D requires 
non-zero values for all bi and all q. Using this, it is straightforward to find the most general 
solution to the set of i^-term equations. Up to symmetry transformations, we obtain 
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(7.10) 



We remark that the trivial vacuum, that is the vacuum configuration where none of 
the flavons develops a VEV, typically provides a solution to the F-term equations as well. 
This can be eliminated by including soft supersymmetry breaking effects. Then the scalar 
potential relevant for the flavon alignments takes the general form 
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(7.11) 



where m^. and m^. denote the soft breaking masses of the driving fields Xi and the 
flavons (pi. The dots stand for additional soft breaking terms. Assuming positive m"^,, 
the driving fields do not develop a VEV. As a consequence, the operators which involve 
a driving field, i.e. those represented by the second term of Eq. ( [7.11| ), vanish. The first 
term, on the other hand, only depends on the flavon fields. This together with negative 



removes the trivial vacuum configuration and the flavons acquire a VEV [30,99]. 



Alternatively, it is in principle also possible to add an explicit mass scale in the flavon 
potential which will then drive the flavon VEVs to non-zero values. For instance, the cube 
of the 5*4 doublet flavon 02 of Eq. can be contracted to an ^4 singlet with a non- 

vanishing VEV. Introducing a driving fleld X'^ one could therefore write down the driving 
terms 



X[ 



M 



(7.12) 



where M denotes an explicit mass scale. Since M is a pure (dimensionful) number, the 
driving fleld and with it {4'2)^ must be completely neutral under any imposed symmetry. 
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In particular, they must not carry charges under extra so-cahed shaping symmetries which 
are typicahy introduced in concrete models to separate the flavons of different sectors. 
In the given example, a shaping symmetry under which the flavon 02 (as well as the 
driving field X-j_) carries charge +1 allows for the coexistence of the alignment term of 
Eq. ( |7.5| ) together with the term of Eq. (7.12) which explictly drives the flavon VEV to 
non-zero values. Assuming a CP conserving high energy theory where all parameters of 
the model can be chosen to be real, driving terms of the form of Eq. ( [7. 121 ) could generate 
spontaneous CP violation where the values of the CP violating phases are constrained to 
a finite number of choices [103]. 



7.3 Direct models after Daya Bay and RENO 

The method of identifying the flavon alignments of direct models using Eq. ( [7.1| ) can be 
applied to any mixing pattern. Yet, until recently, the main focus was limited to only a few 
simple cases, namely tri-bimaximal, bimaximal and golden ratio mixing, see Subsection |3.l| , 
all of which predict ^13 = 0° . The observation of a sizable reactor neutrino mixing angle of 
about 9° by the Daya Bay and RENO collaborations in early 2012, preceded by first hints 
for a non-zero ^13 from the T2K, MINOS and Double Chooz experiments in 2011, has now 
ruled out these simple mixing patterns. This fact seems to call into question the direct 
model building approach. However, it is worth recalling that a vanishing reactor angle has 
long been compatible with experimental data, and hence there was no need to consider 
more complicated mixing patterns. The situation has now changed, and new strategies for 
constructing direct models have to be conceived. 

There are two main paths one can pursue. The first is based on small groups like e.g. 
S4 and A^, and leads to the simple TB, BM or GR mixing patterns with vanishing reactor 
angle at leading order. To render such models compatible with a sizable 613 it is critical 
to discuss higher order corrections which break the simple structure of the mixing matrix. 



This situation is depicted in the central branch of Fig. 10 for the family symmetry 84^. 
Depending on which group elements are selected for the symmetries of the charged lepton 
sector (T) and the neutrino sector {S, U), it is possible to obtain either TB or BM mixing 
from 5*4 at leading order. Analogously, the leading order GR mixing derived from 



can be perturbed by higher order effects (not shown explicitly in Fig. IC). In general, 
higher order corrections are guaranteed to perturb the leading order structure by only 
small contributions. The breaking of the leading order structure can happen either in the 
charged lepton or the neutrino sector. The former entails charged lepton corrections of the 
simple leading order mixing patterns, which give rise to solar mixing sum rules as discussed 



in Subsection 3.5. If the breaking occurs in the neutrino sector, it is possible to break either 
one or both Z2 factors of the leading order Klein symmetry. As the U symmetry typically 
enforces 613 = in these models, it is necessary to break U in either case. Demanding S 
to remain a good symmetry at higher order, gives rise to atmospheric mixing sum rules, 



see Subsection 3.6, while breaking also S leads to arbitrary and unpredictive higher order 



corrections. In Subsection 10.2 we will present a concrete ^4 x SU{5) model of tri-bimaximal 



mixing at leading order, augmented by higher order corrections which break U but not S. 
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Figure 10: Possible strategies for constructing direct models after Daya Bay and RENO. Adopting 
small family symmetries G which predict simple LO mixing patters with 6*13 = (e.g. S^^ A^), 
requires HO corrections. Larger family symmetries can give rise to richer LO mixing patterns with 
non-zero ^13 (e.g. A(96)). The A4 family symmetry refers to the semi-direct case as discussed in 
the text. 



This model yields the trimaximal neutrino mixing pattern TM2, see Eq. (3.32), which can 
accommodate a sizable reactor angle. 

The second strategy of constructing direct models compatible with a sizable reactor 
angle makes use of larger groups such as A(96), see left branch of Fig. Such groups are 
capable of predicting richer leading order mixing patterns (e.g. bi-trimaximal mixing [31]) 
as they contain non-standard Klein symmetries, generated by more complicated forms 
of the elements S and/or U [104, 105]. As before, higher order effects can correct these 
leading order mixing patterns. Charged lepton corrections induced by a breaking of the T 
generator give rise to new solar sum rules. Indeed, in the A(96) x SU (5) model discussed in 
Subsection 10.3, the charged lepton corrections are essential in driving the resulting reactor 
angle to a physically viable value. In the neutrino sector, it is generally possible to break 
either S or U , however, in typical models, the symmetry associated with S stabilises the 
solar angle at a phenomenologically viable value. In practice, it should therefore be the U 
generator which gets broken at higher order, leading to new atmospheric sum rules. 

Before turning to the breaking of the family symmetry in indirect models, we comment 
on the case of semi-direct models. As mentioned at the end of Subsection 3.3, the Altarelli- 
Feruglio model [30, 99] provides an example of a semi-direct model. While the tri- 
bimaximal S symmetry of the neutrino mass matrix forms part of the family symmetry, 
the tri-bimaximal U symmetry arises accidentally due to the absence of flavons in the 1' 
and 1" representations of A^. Introducing such neutrino type flavons in the non-trivial one- 
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dimensional representations, a situation which is tantamount to breaking the U generator 
in 5*4 (see Fig. p!o| ), generates contributions to the neutrino mass matrix which are not of 
tri-bimaximal form [102,106]. However, as can be seen from Appendix 1^, they respect the 
original S symmetry, thus enforcing the trimaximal mixing pattern TM2, see Eq. ( |3l3^ ). 
The accidental tri-bimaximal U symmetry, on the other hand, gets broken by the VEVs of 
the flavons in the 1' and 1" representations, which allows to accommodate arbitrary values 
of ^13. The relative smallness of the reactor angle, compared to the solar and atmospheric 
angles, remains unaccounted for and must therefore be understood as a result of a mild 
tuning of parameters. 

A similar semi-direct approach was taken by Hernandez and Smirnov [107] in an ef- 
fort to accommodate a sizable reactor angle. Focusing on the relevant von Dyck groups 
A4, S4 and As, they demand the T symmetry of the charged leptons and (only) one Z2 
factor of the neutrino Klein symmetry to arise as unbroken subgroups of the underlying 
family symmetry. This strategy allows to identify viable mixing patterns in which a given 
column of the PMNS matrix is completely determined by the properties of the imposed 
symmetry group. For the successful cases, these columns are identical (in some cases up 
to permutations of the rows) to either the first or the second column of the bimaximal, 
the tri-bimaximal and the golden ratio mixing patterns, cf. Subsection |3.1| , [108]. With the 
other two columns of the mixing matrix unspecified, the reactor angle can be regarded as a 
free parameter which, together with the CP phase 6, gives rise to predictions for the other 
two mixing angles, expressed in the form of (exact) sum rules. 



8. Indirect model building 

8.1 Flavon alignments in indirect models 

The vast majority of indirect models is formulated in the framework of the type I seesaw 
mechanism where the right-handed neutrino and the charged lepton mass matrices are 



both diagonal, see Subsection |6.4| . The lepton mixing arises from the structure of the 
Dirac neutrino Yukawa matrix, which in turn originates from the alignment of the flavon 
fields (/)•'. With the lepton doublet L furnishing a triplet representation 3 of the family 
symmetry G, the neutrino flavons typically transform as a 3 of G.^^ The family indices are 
then contracted to the G singlet in the familiar SU{3) way, showing that the columns of the 
Dirac neutrino Yukawa matrix are proportional to the alignments of the flavon fields (f)'^, 
as presented in Eq. (|6.14|) . Application of the seesaw formula gives rise to an effective light 



neutrino mass matrix of the form 



m 



LL 



Y,m'^^,<^>f , (8.1) 



=1 



^^If the triplet 3 is real, L and 4)^ transform in the same representation. In indirect models, the basis of 
the triplet representation of G must then be chosen explicitly real. Note that for this reason the Clebsch- 
Gordan coefficients of S4, and A4, given in Appendix ^ are not applicable in indirect models. For a basis 



suitable for indirect models, see Eq. (5.8) and the discussion thereafter 
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where $j oc denotes a dimensionless vector normalised to one. From the model building 
perspective, the direction of these vectors in flavour space depends on the alignment of the 
flavon VEVs. In general one can distinguish two cases. In this subsection we focus on the 
case where the flavon alignments are orthogonal to each other. The situation where this is 
not the case will be treated in Subsection |8.3| . 

Under the assumption that and $j are orthogonal for i ^ j, the light neutrino mass 
matrix m'^]^ of Eq. ( |8.1| ) is diagonalised by a unitary PMNS mixing matrix with columns 

The resulting eigenvalues are simply m?. This scenario, in which the columns of the 
Dirac neutrino Yukawa matrix are proportional to the columns of the PMNS mixing matrix, 
is called form dominance (FD) [100]. An example of this is provided by the alignments of 
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generating the famous tri-bimaximal mixing pattern. It is important to notice that FD 
is a general concept which applies to arbitrary orthogonal vectors <I>j. In principle, one 
could therefore choose the $i directions so as to yield the experimentally observed mixing 
matrix. However, when it comes to building a model of flavour, a crucial ingredient is the 



justification of the assumed flavon alignments, see Subsection B.2. Therefore, in practice. 



only "simple" alignments are adopted in explicit models. With the parameters of 
Eq. ( |8.1| ) being completely independent of the vectors (which arise from some flavon 
alignment mechanism), it is clear that the mixing matrix does not depend on the masses. 
In indirect models, FD thus implies form diagonalisability [100,109].^^ 

A special case of FD is obtained if the three contributions to the neutrino mass matrix 
of Eq. ( ^.l| ) feature a hierarchy m'j* ^ ^ m!]. In such a scenario, which is called 



sequential dominance [51,52,76], see Subsection \i.3[ the first term, and with it the vector 
$1, can be ignored to good approximation. In fact, one can even remove the flavon (p'( 
from the theory altogether. This would set automatically to zero, without affecting the 
pattern of the 3x3 mixing. The latter can be understood by realising that the flrst column 
of the mixing matrix is uniquely determined by requiring orthogonality to the other two 
columns $2 and $3. As above, SD is a general concept applicable to arbitrary two (or 
three) orthogonal flavon alignments. Choosing ^2 and $3 as given in Eq. ( ^.2| ) leads to 
constrained sequential dominance (CSD) [70], and predicts tri-bimaximal neutrino mixing. 



8.2 Vacuum alignment mechanism in indirect models 



We have discussed in Subsection |7.2| how flavons of direct models can be aligned using the 
i^-term alignment mechanism. In indirect models, the same mechanism is available, how- 
ever, if a triplet representation of the family symmetry is real, it is mandatory to work in a 



^'^However, this is not necessarily the case in direct models, see [102], where the columns of the Dirac 
neutrino Yukawa matrix - in the basis of diagonal right-handed neutrinos - are not related to the flavon 
alignments in a simple way. 
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basis where this is exphcitly reahsed, i.e. where ah group generators are real. Applications 
of the F-term alignment mechanism in indirect models can be found e.g. in [103,110,111]. 
In addition to the usual i^'-term alignment mechanism, indirect models offer an elegant al- 
ternative possibility for achieving particular flavon vacuum configurations. This so-called 
D-term alignment mechanism, as the name suggests, was first implemented in supersym- 
metric models [86, 112], however it is also possible to apply it in a non-supersymmetric 
context. 

The starting point is a flavon scalar potential field which may or may not arise in a 
supersymmetric model from L>-terms, 
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+ Ay , 



(8.3) 



where the index i labels the components of a particular flavon triplet cj) and 
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Ignoring the term Ay in Eq. (8.3), the potential features an SU{3) symmetry and, as 
a consequence, no direction of the flavon alignment would be preferred. Inclusion of the 
term Ay breaks the SU{3) symmetry of the potential and leads to minima which single 
out particular vacuum alignments. With the scale of the flavon VEV depending on m^, A 
and K, it is sufficient to consider the extrema of the quartic term in Eq. ( |8.4D for a unit 
vector <I>. If K > 0, it is necessary to minimise the sum I'l'*!^; leading to the solution 
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(8.5) 



where are arbitrary phases. Such an alignment is of the form of ^2 in Eq. ( |8.2| ). In 
fact, in indirect models, where the alignment of Eq. ( ^.5| ) appears as a column of the Dirac 
neutrino Yukawa matrix, the phases -di can be removed by a field redefinition of the charged 
leptons. In the case where k < 0, the sum |$*|^ has to be maximised. This gives rise to 
the alignment 



K < 
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5.6) 



and permutations thereof. Such alignments are typically useful for constructing a diagonal 
charged lepton sector. They are furthermore necessary to obtain the alignments $3 in 
Eq. (S^) via SU{3) invariant orthogonality conditions. Introducing a new fiavon field (j) 
which couples to the flavons (with alignment <I>-|- = $2) and 0_ (with alignment <!>_) 
as 
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(8.7) 
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we generate the alignment {(j)) oc <I>3 if k! and k" are taken to be positive. An alignment 
proportional to of Eq. ( ^.21) can be obtained subsequently from orthogonality conditions 
involving flavons with alignments along the directions ^>2 and $3. 

The preceding discussion illustrates the importance of the SU (3) breaking term in 
Eq. (|8.3| ). It is therefore natural to identify finite groups G which have the operator in 
Eq. ( ^.41 ) as an invariant. Obviously, the family symmetry G must admit at least one 
triplet representation, with generators which are symmetry transformations of Eq. (^.4|). 
As was shown in [97], possible candidate symmetries include the groups A(3n2) [82,85], 
A(6n2) [82,83] and T„ [90], cf. also Eqs. ([5^81,151^). 



All these symmetries allow for at least two quartic invariants of type 3333, namely 
the S'C/(3) invariant and the operator of Eq. ( |8.4| ). However, four of them have additional 
independent quartic invariants. These are A(24) = 5*4 with one extra invariant, as well 
as A(12) = A4, A(27) and A(54) with two additional invariants each [97]. These new 
invariants may spoil the structure of the vacuum derived from AF of Eq. ( ^.4| ) unless they 
are sufficiently suppressed. From this perspective, the groups A(3n^) and A(6n^) with 
n > 3, as well as the groups T„ are preferred candidates for the underlying discrete family 
symmetry of indirect models. 

We conclude the discussion of the alignments in indirect models with a possible alter- 
native to the invariant of Eq. (|8.4| ), which has not received any attention yet. A cubic term 
of the form is left invariant under the groups A(3n^) and T„. As such a term is 

generally not real, the new term in the flavon potential reads 

Ay = + H.c.) , (8.8) 



replacing the operator in Eq. (|8.4| ). One can easily show that, with k < 0, such a term 
in the flavon potential would generate an alignment of type see Eq. ( |8.5| ), with = 

-{'^1+^2). 

8.3 Indirect models after Daya Bay and RENO 

So far, we have only discussed the flavon alignments of indirect models leading to tri- 
bimaximal mixing. The measurement of a reactor angle 613 of around 9° by the Daya Bay 
and RENO collaborations has now ruled out models of accurate tri-bimaximal mixing. 
This fact demands a modification or extension of the common strategies for constructing 
indirect models. 

As for direct models, there are two principle paths one can pursue. The first builds 
on existing indirect models of tri-bimaximal mixing which arise in the framework of con- 



strained sequential dominance, see the central branch of Fig. 11. Such a leading order 
structure must be broken by higher order corrections, which can stem from either the 
charged lepton or the neutrino sector. The former case requires a breaking of the acci- 
dental T symmetry and leads to solar mixing sum rules as discussed in Subsection |3.5| . 
Alternatively, the higher order corrections can break the accidental tri-bimaximal U sym- 
metry in the neutrino sector, entailing atmospheric mixing sum rules, see Subsection |3.6| . 
Breaking the tri-bimaximal Klein symmetry of the neutrino sector completely gives rise 
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Figure 11: Possible strategies for constructing indirect models after Daya Bay and RENO. Starting 
from CSD it is possible to add HO corrections. Alternatively, it is possible to modify the flavon 
alignments of CSD at leading order, thus giving rise to CSD2 and PCSD. 

to arbitrary and therefore unpredictive corrections to the mixing angles. An example for 
higher order corrections which break the U symmetry can be constructed using the flavon 
alignments of constrained sequential dominance, proportional to $2 and $3, see Eq. (^.2|), 
and add a small perturbation along the direction $1 to $3, see also [113]. The form of the 
two flavon alignments can then be written as 

^ oc $2, ^ oc c|,'3 = $3 + ecI>i, (8.9) 

with e ^ 1. Note that these two vectors are still orthogonal to each other, hence the 
conditions of form dominance are satisfied. With these alignments, it is straightforward to 
show that the resulting light neutrino mass matrix 

<L = "Z^^2^^ + "^3^3^f , (8.10) 

breaks the original U symmetry while continuing to respect the S symmetry. The latter 
can be easily verified by noticing that (1, 1, 1)"^ is still an eigenvector of m'^^ in Eq. ( 8.10| ), 



meaning that the second column of the tri-bimaximal PMNS mixing matrix, i.e. ^2, re- 
mains unchanged. Hence the trimaximal mixing structure TM2, see Eq. ( p. 32 ), is achieved, 



which stabilises the solar angle. On the other hand, the breaking of the U symmetry of 
the neutrino mass matrix by higher order effects allows to accommodate non-zero ^13. 

The second strategy of constructing indirect models with sizable ^13 is based on new 
alignments at leading order. In the following we present two examples (see the right and 
the left branch of Fig. |l^): partially constrained sequential dominance (PCSD) [73] and 
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constrained sequential dominance 2 (CSD2) [111]. Both scenarios make use of two flavon 
triplets whose alignments are not orthogonal to each other, in contrast to the previously 
discussed cases. 



PCSD 



Partially constrained sequential dominance was first proposed in [73] as a simple modifica- 
tion of CSD, where one flavon is aligned along the original $2 direction, while the alignment 
of the other flavon is assumed to deviate slightly from the $3 direction by filling the zero 
of the first component with e <C 1, i.e. 
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Note that $2 and <1>3 are not orthogonal to each other, hence PCSD violates form dominance 
at linear order in e. Inserting these two alignments into Eq. (|8.lD yields the effective 
neutrino mass matrix 

fill' 
1 1 1 
1 1 1, 
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(8.12) 
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This matrix - with non-zero e - is no longer diagonalised by a the tri-bimaximal mixing 
matrix. Assuming |e| w 0.2 as well as a normal neutrino mass hierarchy, i.e. \m2\ ~ le^^sl; 
analytic expressions for the mixing parameters valid to second order in e were derived 
in [114]. These results show that to first order in e, the tri-bimaximal solar and atmospheric 
mixing angle predictions are maintained while the reactor angle takes a value of order e. 
Therefore, PCSD gives rise to tri-bimaximal-reactor mixing, see Eq. (|]2|), at leading 
order. A special case of TBR mixing is obtained if the parameter e can be identified with 
the Wolfenstein parameter A = 0.2253ib0.0007. As discussed in [67], such a situation results 
in a reactor angle which satisfies sin^is = leading to 6*13 9.2°, a value remarkably 



close to the one measured by Daya Bay and RENO, see Subsection ^ 
The alignment of the fiavon 



63 in Eq. ( |8.11| ) can be achieved through both the F- 
term and the D-term alignment mechanism. The starting point are the simple alignments 
proportional to 



1 

7i 



1 



.13) 



v-v 



as obtained for instance in Eq. ( |8.2| ) and Eq. (|8.6[ ), respectively. The alignment in the 
direction of $3 then arises from successive orthogonality conditions as follows. Imposing 
orthogonality of the VEV of an auxiliary fiavon (j)a with (cp^) and {(px) yields 

/o\ 

± ((J)-^) and ((j)r,) ± U.) ^ « — 1 . (8.14) 
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Requiring the alignment of the flavon (p^ to be orthogonal to the alignment of this auxiliary 
flavon, we find the general structure 



_L 



A 



oc 




(8.15) 



where rii and n2 can take arbitrary values. For ni <^ n2, this is nothing but the alignment 
in the direction of ^3. The hierarchy between ni and n2 may either be a consequence of 
mild tuning or, under certain assumptions, result from a combination of a renormalisable 
and a non-renormalisable term, where, after contracting the family indices, the former is 
proportional to rii while the latter is proportional to n2. The necessary mass suppression 
of the non-renormalisable term then naturally suppresses e = ^ ^ 1 [114]. In order to 
establish a connection of e and A, one can envisage scenarios in which the flavon 03 appears 
in both the neutrino as well as in the quark sector. In the latter, cp^, has to be responsible 
for generating the Cabibbo mixing [67]. 



CSD2 

Constrained sequential dominance 2, proposed in [111], assumes two flavon fields in the 
neutrino sector. One is aligned along the direction of ^3 of Eq. ( ^.21) , while the alignment 
of the other flavon is a relatively simple variation of ^2 of Eq. (|8.2|), explicitly 
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or 
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-1, 



Analogous to the case of PCSD, the alignments of the two flavons, pointing in the directions 
of and $3, are not orthogonal to each other, implying that CSD2 violates form domi- 
nance as well. In the following, we only present the discussion of the first ^"2 alignment of 
Eq. ( ^.16|) ; the alternative case of the second alignment can be treated analogously, leading 
to almost identical results. Inserting the alignments of (j)2 and (ps into Eq. (^?^) generates 
the neutrino mass matrix 



/l2 0) 
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(8.17) 



which is not of tri-bimaximal structure. Yet, one can immediately verify that (—2, 1, 1)"^, 
i.e. the first column of the tri-bimaximal mixing matrix, is still an eigenvector of m'^^. This 
shows that CSD2 necessarily leads to the trimaximal mixing pattern TMl, see Eq. ( 3.30| ). 



With the assumption of a normal neutrino mass hierarchy, i.e. with |m2| ~ k^isl) analytic 
expressions for the mixing parameters valid to second order in e were derived in [111]. 
These results explicitly confirm that the solar angle maintains its tri-bimaximal value at 
linear order in e, while the deviations of the reactor and the atmospheric mixing angles 
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from their tri-bimaximal values are proportional to e, leading to the linear mixing sum rule 
of Eq. ( pl| ). 

The alignment of 02 in Eq. ( |8.16| ) can be derived from a set of orthogonality conditions 
similar to the situation in PCSD. In CSD2 we start from the simple alignments proportional 
to 
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see Eqs. (p^3-€)- Demanding orthogonality of the flavon VEV (^2) with ((/>i) and 
generates an alignment in the direction of the first <I>2 vector in Eq. (p.l6), 
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and 
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The second $2 vector in Eq. (^.16|) can be obtained similarly by using instead of (^z- 



9. Grand unified theories of flavour 



9.1 Grand unified theories 

One of the exciting things about the discovery of neutrino masses and mixing angles is 
that this provides additional information about the flavour problem - the problem of un- 
derstanding the origin of three families of quarks and leptons and their masses and mixing 
angles. In the framework of the seesaw mechanism, new physics beyond the Standard 
Model is required to violate lepton number and generate right-handed neutrino masses 
which may be as large as the GUT scale. This is also exciting since it implies that the 
origin of neutrino masses is also related to some GUT symmetry group Ggut, which unifies 



the fermions within each family as shown in Fig. 12. Some possible candidate unified gauge 



groups are shown in Fig. |13. 

Let us take Gqut = 5'^^ (5) as an example. Each family of quarks (with colour r, h, g) 
and leptons fits nicely into SU{5) representations of left-handed (L) fermions, F = 5 and 
T = 10 

/O < - 



d'^ 

e~ 



n: 







Ug dg 





(9.1) 



0/ 



where c denotes CP conjugated fermions. The SU{5) representations F = 5 and T = 10 
decompose into multiplets of the SM gauge group SU{3)c x SU{2) l^U (l)y as F = (d^, L), 
corresponding to, 

5 = (3,1,1/3)0(1,2,-1/2), (9.2) 
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Figure 12: The fermion masses are here represented by a lego plot. We have multiplied the masses 
of the bottom, charm and tau by 10, the strange and muon by 10^, the up and down by 10^, the 
electron by 10^ to make the lego blocks visible. Assuming a normal neutrino mass hierarchy, we 
have multiplied the third neutrino mass by 10^"'^ and the second neutrino mass by 10^^ to make the 
lego blocks visible. This underlines how incredibly light the neutrinos are. The symmetry groups 
Ggut and Gfam act in the indicated directions. 

and T = {u^, Q, e^), corresponding to, 

10 = (3, 1, -2/3) e (3, 2, 1/6) © (1, 1, 1). (9.3) 

Thus a complete quark and lepton SM family {Q, u^, dP, L, e*^) is accommodated in the 
F = 5 and T = 10 representations, with right-handed neutrinos, whose CP conjugates are 
denoted as u'^, being singlets of SU (5), i^^ = 1. The Higgs doublets and which break 
electroweak symmetry in a two Higgs doublet model are contained in the SU (5) multiplets 
and H-^. 

The Yukawa couplings for one family of quarks and leptons are given by, 

VuH^iT^kTime'^'''^ + yuH^iF'u" + VdHiTi.F^, (9.4) 

where e*-?'^''" is the totally antisymmetric tensor of SU{5) with k,l = 1, . . . , 5, which 
decompose into the SM Yukawa couplings 

VuHuQu^ + y.HuLu^ + yd{HdQ<r + H^e'^L). (9.5) 
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SU{5) X Uil) 




SUi3)c X SUi3)L X SUi3)R 



SUiA)c X SU{2)l X SUi2)R 



SUi3)c X SUi2)L X 5[/(2)r X 



^^/{a)^ X 5C7(2)i X U{1)y 



Figure 13: Some possible candidate unified gauge groups. 

Notice that the Yukawa couphngs for down quarks and charged leptons are equal at the 
GUT scale. Generalising this relation to all three families we find the SU{5) prediction for 
Yukawa matrices, 

Yd = Y^, (9.6) 

which is successful for the third family, but fails badly for the first and second families. 
Georgi and Jarlskog [115] suggested to include a higher Higgs representation H-^ which is 
responsible for the 2-2 entry of the down and charged lepton Yukawa matrices. Dropping 
SU{5) indices for clarity, 

(1^)22^45^2^2, (9.7) 

decomposes into the second family SM Yukawa couplings 

{Yd)22{HdQ2d'i - 3Fde^L2), (9.8) 

where the factor of —3 is an SU{5) Clebsch-Gordan coefficient.^^ Assuming a hierarchical 
Yukawa matrix with a zero Yukawa element (texture) in the 1-1 position, results in the 
GUT scale Yukawa relations, 

Vb = Vt, Vs = y , yd = 3ye, (9.9) 
^^In this setup, Hd is the light linear combination of the electroweak doublets contained in H-^ and -H45-. 
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which, after renormahsation group running effects are taken into account, are consistent 
with the low energy masses. The precise viabihty of these relations has been widely dis- 
cussed in the light of recent progress in lattice theory which enable more precise values 
of quark masses to be determined, especially the strange quark mass (see, e.g., [116]). In 
super symmetric (SUSY) theories with low values of the ratio of Higgs vacuum expectation 
values, the relation for the third generation yi, = Ur at the GUT scale remains viable, 
but a viable GUT scale ratio of Vfi/ys is more accurately achieved within SUSY SU{5) 
GUTs using a Clebsch factor of 9/2, as proposed in [117], which is 50% higher than the 
Georgi-Jarlskog prediction of 3. 

9.2 Combining GUTs and family symmetry 

As already remarked in Subsection p.3| , it is a remarkable fact that the smallest leptonic 
mixing angle, the reactor angle, is of a similar magnitude to the largest quark mixing 
angle, the Cabibbo angle, indeed they may even be equal to each other up to a factor of 
\/2. Such relationships may be a hint of a connection between leptonic mixing and quark 
mixing, where such a connection might be achieved using GUTs [118,119]. For example, the 
Georgi-Jarlskog relations discussed above already lead to the left-handed charged lepton 
mixing angle having a simple relation with the right-handed down-type quark mixing angle 

^ ^12/3 where the approximation assumes hierarchical Yukawa matrices, with the 1-1 
elements being approximately zero. If the upper 2x2 Yukawa matrices are symmetric (as 
motivated by the successful Gatto-Sartori-Tonin (GST) relation [120] which relates the 12 
mixing 9^2^ to the down and strange mass by ^^2'^ ~ y/^fn^Jm^) then we may drop the 
L,R subscripts and this relation simply becomes 9^2 = ^i2/3- In large classes of models, 
including those discussed later, the quark mixing originates predominantly from the down- 
type quark sector, in which case this relation becomes 9^2 = If one starts from TB 
mixing in the neutrino sector, resulting from some discrete family symmetry, then, using 
the results in Subsection |3.4| such a charged lepton correction results in a reactor angle in 
the lepton sector of ^13 ~ 9c/{3\/2) as discussed for example in [70]. This is a factor of 3 
too small to account for the observed reactor angle, but it illustrates how the reactor angle 
could possibly be related to the Cabibbo angle using GUTs. Indeed it has been suggested 
that perhaps the charged lepton mixing angle is exactly equal to the Cabibbo angle in some 
GUT model, leading to 6*13 ^ 9c /V^ [67, 72, 121]. However it is non-trivial to reconcile 
such large charged lepton mixing with the successful relationships between charged lepton 
and down-type quark masses, and it seems more likely that charged lepton mixing is not 
entirely responsible for the reactor angle. 

The above discussion provides an additional motivation for combining GUTs with 
discrete family symmetry in order to account for the reactor angle. Putting these two 
ideas together we are suggestively led to a framework of new physics beyond the Standard 
Model based on commuting GUT and family symmetry groups, 

Ggut X Gfam. (9.10) 

There are many possible candidate GUT and family symmetry groups some of which are 
listed in Table |2[ Unfortunately the model dependence does not end there, since the 
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Table 2: Some candidate GUT and discrete family symmetry groups, and the papers that use 
these symmetries to successfully describe the solar, atmospheric and reactor neutrino data. 

details of the symmetry breaking vacuum plays a crucial role in specifying the model and 
determining the masses and mixing angles, resulting in many models. These models may 
be classified according to the particular GUT and family symmetry they assume as shown 
in Table |[ 

Another complication is that the masses and mixing angles determined in some high 
energy theory must be run down to low energies using the renormalisation group equa- 
tions. Large radiative corrections are seen when the seesaw parameters are tuned, since 
the spectrum is sensitive to small changes in the parameters, and this effect is sometimes 
used to magnify small mixing angles into large ones. 

In natural models with a normal mass hierarchy based on SD the parameters are 
not tuned, since the hierarchy and large atmospheric and solar angles arise naturally as 
discussed in the previous section. Therefore in SD models the radiative corrections to 
neutrino masses and mixing angles are only expected to be a few per cent, and this has 
been verified numerically. 

10. Model examples 

In this section we give three examples of GUTs of flavour based on the (semi-)direct 
approach, which can account for all quark and lepton masses and mixing, including the 
observed reactor angle which only gets a small correction from charged lepton mixing. The 
first example is based on the minimal family symmetry A4 combined with the minimal 
GUT SU{5). This is actually a semi-direct model, since only half the Klein symmetry is 
contained in A4, resulting in TM2 mixing, but, like all semi-direct models, it cannot explain 
the relative smallness of the reactor angle. The second example does explain the smallness 
of the reactor angle compared to the atmospheric or solar angles by embedding the A4 into 
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Table 3: The charge assignments of the matter, Higgs and flavon superfields in the A\ x SU{5) 
model of [132]. The shaping symmetry U{1) x Z2 x Z3 x Z5 constrains the set of operators allowed 
in the superpotential. 



an S4 family symmetry. This allows a direct model where the Klein symmetry resulting 
from S4 is half broken by a rather special higher order correction, resulting again in TM2 
mixing as in the case. The third model based on A (96) is an example of a direct model 
with a larger family symmetry where the Klein symmetry corresponds to bi-trimaximal 
mixing in which the reactor angle is already non-zero at the leading order, and where the 
small charged lepton correction with an assumed zero phase brings it into agreement with 
the Daya Bay and RENO measurements of ^13 ~ 9°. 

10.1 A4 X SU{5): a semi-direct model of trimaximal mixing 

Our first example of a Gfam x Ggut model with large ^13 is based on the Altarelli- 
Feruglio model of leptons [30,99]. Working in a super symmetric SU{5) setting, the 



three matter families of -F = 5 and T = 10, see Subsection |9.1| , transform under A4 
as 3 and l",l',l, respectively. The seesaw mechanism is implemented in the model by 
introducing right-handed neutrinos v'^ living in the 3 of A4. The Higgs fields H^, H-^ and 
^45 furnish the one-dimensional A4 representations 1, 1' and 1". The latter gives rise to 
the intriguing Georgi-Jarlskog (GJ) relations [115]. The family symmetry breaking flavon 
fields are SU (5) singlets and can be divided into fields which appear in the neutrino sector 
(p'^ and fields which appear in the quark sector (p^. The ^4 family symmetry is enriched by 
the shaping symmetry f7(l) x Z2 x Z3 x Z5 in order to control the coupling of the flavons 
to the different matter sectors. The complete charge assignments of the matter, Higgs and 
flavon superfields is presented in Table ^. 

A discussion of all the different aspects of the model, including the vacuum alignment, 
can be found in [132]. Here we mainly focus our attention on the neutrino sector. The 
leading order operators of the neutrino superpotential which are invariant under all imposed 
symmetries of the model, see Table ^ reads 

~ Fi^'^Hr, + (0g + (1)1 + <t>l, + ^l„) v^v^ , (10.1) 
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where dimensionless order one coupling coefficients are suppressed. Note that the Dirac 
Yukawa term does not involve any ffavon field, hence, the Dirac neutrino Yukawa matrix 
does not break the symmetry. Inserting the Higgs VEV Vu (the VEV of the electroweak 
doublet contained in H^) then generates an A4 invariant Dirac neutrino mass matrix rriLfi. 
As a consequence, the mixing in the neutrino sector originates solely from the right-handed 
neutrino mass matrix whose form depends on the chosen basis of as well as the vacuum 
configuration of the flavon fields. Working in the A4 basis of Appendix j^, it has been 
shown in [132] that the two triplet flavon alignments 




(10.2) 



can be obtained from the i^-term alignment mechanism along the lines of [30]. This ne- 



cessitates the introduction of a t/(l)_R symmetry, a set of driving fields, see Subsection 7^ 
as well as the auxiliary flavon field (pl which has the same charges as the triplet flavon 
(/)3, except for A4. With these alignments as well as VEVs for the one-dimensional flavon 
fields, the Dirac neutrino mass matrix and the right-handed neutrino mass matrix take the 
form^^ 

/l o\ 

niLR ^ 1 , (10.3) 
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The effective light neutrino mass matrix emerges from these by applying the type I seesaw 
formula of Eq. (1^), i.e. 

"^LL = -mLRM]^^mlji . (10.5) 

Due to the trivial structure of the Dirac neutrino mass matrix rriLji, the neutrino mixing 
matrix is identical to the unitary matrix which diagonalises Mrji (and automatically also 
M^^), except for a permutation of the second and the third row [102]. We remark that the 
light neutrino masses are not related by a mass sum rule since the right-handed neutrino 
mass term involves four independent flavon fields. The inverse mass sum rule for A4 quoted 
in Eq. (^]^) with 7=1 and 6 = —2 can only be recovered by removing the flavons (f)'^, and 

, which in turn corresponds to the well-known case of tri-bimaximal neutrino mixing. 

Rather than diagonalising Mrr explicitly, let us discuss its symmetries. The first two 



terms in Eq. ( 10.4 ) are symmetric under the tri-bimaximal Klein generators S and U of 



^As the mass term of the right-handed neutrinos in the superpotential of Eq. (10.1) involves the CP 



conjugate fields u"^, the mass matrix in the conventions of Subsection 4.1 is obtained using the complex 



conjugate flavons VEVs ^Pr* , cf. also Eq. (i.l3). However, for notational clarity, we drop the * here and in 
the following. 
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Eq. (|6.^). The third and the fourth terms break the tri-bimaximal structure, however, 
in a special way. It is straightforward to prove exphcitly that S'^ MrrS = Mrr is stiU 
respected. A simple way of seeing this is by noticing that all neutrino flavon VEVs remain 
unchanged under the ^4 transformation S, see Appendix |y. On the other hand, the U 



matrix of Eq. (3.8) does not form part of A4. In order for M/j^ to be also symmetric under 
[/, hence entailing tri-bimaximal neutrino mixing, one would have to require (py = fyi- 
In [30, 99] this condition among the a priori unrelated VEVs is realised by not including 
flavons in the 1' and 1" representations of A/^ in the first place. Alternatively, the two non- 
trivial one-dimensional representations can be unified into a doublet of 5*4, see Appendix^. 
A suitable VEV alignment of such a doublet can relate the two components such as to 
generate a right-handed neutrino mass matrix of tri-bimaximal form, see e.g. [102]. In 
general, however, ^ ip^,, and there is no accidental U symmetry in an ^4 model with 
neutrino flavons in all possible representations of the family symmetry. 

With nriLR being invariant under the full ^44 family symmetry and Man being symmet- 
ric under S*, the type I seesaw mechanism generates a light effective neutrino mass matrix 
which also respects the S symmetry. This symmetry can be translated to a particular 
mixing pattern by considering an eigenvector v oi S with eigenvalue +1. One can easily 
check that the only solution to Sv = 7; is of the form v oc (1, 1, 1)^. Using this and the 
invar iance of the light neutrino mass matrix m-^j^ under S as well as S'^ = S, we obtain 

^IlV = StullSv = Sml^v , (10.6) 

which shows that m'^j^v is an eigenvector of S with eigenvector +1, and thus 

m'^i^v oc V . (10-7) 

As V oc (1, 1, 1)"^ is an eigenvector of the neutrino mass matrix, the normalised vector ^ 
corresponds to a column of the neutrino mixing matrix. Except for being orthogonal to v, 
the other two columns are not specified by the S symmetry. In order to be meaningful for 
physics, the vector has to be identified with the second column of the neutrino mixing 
matrix, so that the trimaximal pattern TM2 of Eq. ( |3.32| ) ensues. As mentioned earlier, 
this special mixing pattern allows for a large reactor angle while keeping the solar angle at 
its tri-bimaximal value at leading order. 

It is now important to notice that the TM2 mixing sum rule of Eq. ( p. 33 ) applies only 



to the neutrino sector. In other words, the TB deviation parameters as well as the CP 
phase should carry a neutrino index, i.e. s'^,a'^,r'^ and S'^ . In order to find the deviation 
parameters of the physical PMNS matrix, it is necessary to add the effect of the charged 
lepton corrections. As is demonstrated in [132], the mass matrices of the charged fermions 
in the A4 x SU{5) model are given by 



/e6 g4 



e 



/-3e6 

-3e^ -3e'' I v, , (10. 



V 




where the scales of the flavon VEVs were assumed to be 

(/J? ~ ~ eM , ~ e^M . (10.9) 
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M denotes a messenger mass which we allow to vary for the up-type and the down-type 
quarks, thus justifying a different expansion parameter in (e) and M^^ e (e)- The factors 
of —3 in Me correspond to an SU{5) Clebsch-Gordan coefficient which originates from 
Georgi-Jarlskog terms of the form FTH-^ [115], multiplied by appropriate products of 
flavon fields. Finally, Vd denotes the VEV of the light combination of the electroweak 
doublets contained in H-^ and H^. 

With the structure of Mg given in Eq. ( |10.^ ), the only significant left-handed charged 
lepton mixing is the 12 mixing ~ e/3. The parameter e can be approximated by 
the Wolfenstein parameter A as the effect of the left-handed up-type quark mixing on the 
CKM matrix is negligible. Combining the TM2 mixing of the neutrino sector and charged 



lepton corrections with ~ A/3, see Subsection 3.4, leads to the sum rule bounds [132] 



|a|<^(r + ^)|cos<5| , \s\<^, (10.10) 

where s, a, r are the physical TB deviation parameters of the PMNS matrix and 5 denotes 
the physical CP phase. 

We conclude the discussion of the x SU (5) model of trimaximal neutrino mixing 
by pointing out that this framework does not provide any explanation for the suppression 
of the reactor angle compared to the solar or atmospheric angles. Therefore, this model 
relies on mild tuning of parameters. In the next subsection we show how to obtain such a 
suppression in the context of an ^4 model of tri-bimaximal mixing in which the U symmetry 
gets broken by higher order corrections. 



10.2 5^4 X SU{^): a direct model of tri-bimaximal mixing vi^ith corrections 

In this subsection we present the main ingredients of the supersymmetric 54 x SU{5) 
model of [134]. It is based on an earlier direct model [137] which has been ruled out by the 
measurement of ^13 w 9°. In order to accommodate this experimental result, the model 
of [137] has simply been augmented with an extra 5*4 singlet flavon field rj. The three 
families of SU{5) matter multiplets -F = 5 and T = 10 transform under 6*4 as 3 and 2-1-1, 
respectively. We furthermore introduce three right-handed neutrinos u'^ which are unified 
in the 3 of S4 and allow for the type I seesaw mechanism. The Higgs sector is 5*4 blind and 
comprises the standard SU{5) Higgses in the 5 and 5, plus an additional Georgi-Jarlskog 
Higgs in the 45. The family symmetry is broken by a set of flavon fields transforming 
in various representations of 5*4. In order to control the coupling of the flavon fields to 
different matter sectors, we impose a global U{1) shaping symmetry. The complete charge 
assignments of matter, Higgs and flavon fields are listed in Table ^ 

With the model formulated at the effective level, it is straightforward to derive the 
leading operators of the matter superpotential which are invariant under all imposed sym- 
metries. Assuming a generic messenger mass M of order the GUT scale, and suppressing 
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Table 4: The charge assignments of the matter, Higgs and flavon superfields in the S4 x SU{5) model 
of [134]. The U{1) shaping symmetry constrains the set of operators allowed in the superpotcntial. 



all dimensionless order one coupling coefficients, we find 

W ~ TsT3H^ + ^TTcl>^H^ + ^TT<P^^^H^ (10.11) 
+^FT3</.i%+ ^(F^i)i(r</.^)ii/35 + ^(F<A^0^)3(r^i)3% (10.12) 



M 



'4v 



(10.13) 



The last term of Eq. ( |10.13 ), highlighted in red colour, provides the source of the higher 
order correction to the right-handed neutrino mass matrix which is essential in generating 
a large reactor angle. In principle, all independent invariant products of the ^4 represen- 
tations have to be considered for each of these terms; in practice, there is often only one 
possible choice. In our example, the second and the third term of Eq. ( 10.12[ ) would give 
rise to several independent terms. However, the contractions specified by the subscripts 
1 and 3 single out a unique choice. Within a given UV completion, the existence and 
non-existence of certain messenger fields can justify such a construction. 

The Yukawa matrices are generated when the flavon fields acquire their VEVs. The 
explicit form of these matrices depends on the S4 basis which we choose as given in Ap- 
pendix]^. Adopting the F-term alignment mechanism which requires to introduce a U{l)ji 



symmetry as well as new driving fields, see Subsection 7.2, is has been shown in [134, 137] 
that the following alignments can be obtained, 



(4) 





^1- 



(10.14) 
(10.15) 

(10.16) 



Inserting these vacuum alignments and the Higgs VEVs Vu and Vd yields a diagonal up- 
type quark mass matrix ~ diag ((/52^2/-^^; ¥^2/-^ 5 1) ^« well as down-type quark 
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and charged lepton mass matrices 
/ 








{^2 



Vd 







Vd ■ 



(10.17) 



(10.18) 



The factors of —3 in Mg originate from the second term of Eq. ( |10.12| ) involving the Georgi- 
Jarlskog Higgs field H-^ [115]. Note that the 1-2 and 2-1 entries, which originate from the 
same superpotential term, have identical absolute values; together with the zero texture 
in the 1-1 entry, this allows for a simple realisation of the GST relation in the 5*4 x SU{5) 
model. In the neutrino sector we find the Dirac neutrino mass matrix and the right-handed 
neutrino mass matrix 

/ 




Mrr 



^1 



- ^3' + 



M 
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(10.19) 

It is clear from Eqs. (10.17- r0.19| ) that the fermion masses and mixings are solely deter- 
mined by the scales of the flavon VEVs. In order to achieve viable GUT scale hierarchies 
of the quark masses and mixing angles [116], we have to assume 

(10.20) 
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where A denotes the Wolfenstein parameter. With these magnitudes, the charged fermion 
mass matrices are fixed completely, 
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(10.21) 



/ A^ A^' 
A^ A^ A^ 
^0 A^ 

Due to the GJ factor of —3 and the texture zero in the 1-1 entry, we obtain viable charged 
lepton masses. With the vanishing off-diagonals in the third column of Mg, there is only 
a non-trivial 12 mixing in the left-handed charged lepton mixing Ve^ , see Subsection |3.4| . 
This mixing, ~ '^/3, will contribute to the total PMNS mixing as a charged lepton 
correction. 

Turning to the neutrino sector, we first observe that the Dirac neutrino Yukawa term 
does not involve any flavon field. As the family symmetry S4 remains unbroken by itilr, 
the mixing pattern of the effective light neutrino mass matrix m'£j^ (obtained from the 
type I seesaw mechanism) is exclusively determined by the structure of Mrr. Dropping 
the higher order terms which are written in red, we note that the leading order structure of 
Mrr, and with it m'^j^, is of tri-bimaximal form.^^ This can be easily seen by verifying that 

^''Similar to the A4, x SU{5) model of Subsection 10.1, the masses of the light neutrinos are not related 
by any mass sum rule as the right-handed neutrino mass matrix AIhr is generated from the VEVs of three 
independent flavon fields. 
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the flavon ahgnments of Eq. ( |10.16| ) are left invariant under the S and U transformations 
of Appendix ^ This leading order tri-bimaximal structure yields light neutrino masses 
m'^ ~ 0.1 eV if we set ^123' ~ A'^M. As we want to break the TB Klein symmetry 
by means of the flavon r] at higher order, we set (77) ~ A^M. Then the TB breaking 
effect is suppressed by one power of A compared to the leading order. The effective flavon 
(j)2'i] transforms as an S4 doublet with an alignment proportional to (1,0)^. As can be 
seen from the S4 generators of the doublet representation, see Appendix Q this alignment 
breaks the U symmetry but respects S. This directly proves that Mjiji as well as m'^j^ 
are both invariant under S, which in turn entails the TM2 neutrino mixing pattern, where 
the second column of the mixing matrix is proportional to (1,1,1)^, cf. Eq. ( |3.32 ). The 
physical PMNS matrix is obtained from multiplying the TM2 neutrino mixing with the 
left-handed charged lepton mixing, see Eq. ( f3.15 ). As a result, we find the same sum rule 
bound as in the A4 x SU{5) model of Subsection 10.1 , given explicitly in Eq. ( |10.10D [134]. 

In summary, the measurement of large ^13 has ruled out the original 5*4 x SU (5) 
model [137] which predicted accurate tri-bimaximal neutrino mixing plus small charged 
lepton corrections. A modest extension of the particle content can induce a breaking of the 
U symmetry of the TB Klein symmetry at relative order A. The required new flavon field 
allows for large 6*13 and does not destroy the successful predictions of the original model, 
i.e. it does not have any significant effects on the quark or flavon sectors of the model. 

10.3 A(96) X SU{5): a direct model of bi-trimaximal mixing 

The direct model discussed in this subsection is based on the observation that larger family 
symmetry groups can contain physically interesting Z2 x Z2 subgroups which, in the basis 
of diagonal charged leptons, differ from the well-known TB Klein symmetry [104,105]. The 
first model of leptons adopting the family symmetry A(96) was constructed in [127]. Here 
we present the first (supersymmetric) grand unified model based on A(96) x SU{5) [31]. 

The group A(96) is a member of the series of groups A(6n^) [83] with n = 4. Like its 
subgroup S4 = A(6 • 2^), it can be obtained from three generators 5, T, U. The group has 
ten irreducible representations: 1, 1', 2, 3, 3, 3, 3', 3', 3', 6. The generators of the one- 
and two-dimensional representations are identical to the corresponding ^4 representations, 
see Appendix |C|. The triplet 3 is generated by 

f-i + Vs-i-Vs -1 

-1-V3 -1 -l + \/3 
-1 -l + V3-l-\/3y 

(10.22) 

where oj = e~ . Notice that S is identical to the tri-bimaximal S symmetry of Eq. (|6. j). 
Invariance of the neutrino sector of a A (96) model under S therefore implies TM2 neutrino 
mixing. The complex conjugate representation 3 is generated by the matrices of Eq. (|10.22| ) 
with T — )• T* . The third triplet 3 is obtained from 

/l 0^ ^ 
Olol, C/ = 4-2 1-2|, T= 010. (10.23) 
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Table 5: The charge assignments of the matter, Higgs and flavon supcrficlds in the A(96) x SU{5) 
model of [31]. The U{1) shaping symmetry is defined by four independent integers x, y, z, and w. 



The generators of the three representations 3, 3, 3 all have determinant +1. This is not 
so for the other three representations 3', 3', 3' which can be obtained from the unprimed 
triplets by simply changing the overall sign of the corresponding U generator. Concern- 
ing the sextet representation and for more details on the group theory of A(96) such as 
Kronecker products and Clebsch-Gordan coefficients, we refer to the extensive appendix 
of [31]. 

The construction of the A (96) x SU (5) model follows closely the logic of the ^4 x SU (5) 
model in [137]. In particular, the flavons of the up-type and the down-type quark sector 
are almost identical, and the GJ mechanism is also implemented along with the GST 
relation. The complete charge assignments of the matter, Higgs and flavon fields of the 
A(96) model [31] are listed in Table ^. In addition to a U{1) shaping symmetry - defined 
in terms of suitably chosen integers x,y,z,w - a Z3 factor has been introduced to forbid 
dangerous terms in the superpotential which, otherwise, would be allowed. 

With the particle content and the symmetries of Table ^, the leading order operators 
of the matter superpotential take the form 



W ~ T^TsH^ + ^TTcb^H^ + ^TTcb^^^H^ (10.24) 
+^Fr34%+ ^(F^^)i(T</.i)i//45 + ^(F0^<j(>^)3(T^)3% (10.25) 
+Fu'H^ + u'^u'^ct)'^, + i/V^^, . (10.26) 



Dimensionless order one coupling coefficients are suppressed, and M is a generic messenger 
mass scale. The subscripts on the parentheses denote the specific contractions being taken 
from the A (96) tensor product contained inside the parentheses. 

As has been elaborated in [31] , the F-term alignment mechanism can be used to derive 
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the vacuum alignment of the flavon fields as 
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(10.27) 
(10.28) 

(10.29) 



Inserting the flavon and the Higgs VEVs yields the charged fermion mass matrices 
diag i^^^^/M^ ip^M , 1) vu and 
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(10.31) 



In the neutrino sector we obtain the Dirac neutrino mass matrix and the right-handed 
neutrino mass matrix 
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(10.32) 

The mass matrices M„, Md and Me of the A(96) model are almost identical to those of the 



54 model presented in Subsection |10.2| . Indeed, with a similar hierarchy of flavon VEVs, 



(^2 ~ 



^1 



(^2 ~ 



AM 



(10.33) 



one obtains the charged fermion mass matrices of Eq. ( 10.21| ). As in Subsection 10.2 
the left-handed charged lepton mixing is described by a unitary matrix with angles 



^23 



and ~ A/3 as well as a general phase SI2, see e.g. Eq. (| 



In the neutrino sector, the right-handed neutrino mass matrix Mjm involves the VEVs 
of only two flavon fields. The three eigenvalues of Mrr are therefore related, leading to the 
mass sum rule for the light neutrinos reported in Eq. (|1.5|) with 7 = 1 and 6 = ±21. Mjm 



^^The alignments presented in Eq. ( 10.29|) do not break the neutrino Klein symmetry generated by S 
and U. While {(p^i) corresponds to the most general such alignment of a flavon in the 3 , this is not the 
case for {<?f>|,). It is straightforward to verify that the most general alignment of a flavon in the 3' of A (96 ) 
which satisfles S'(</)~,) = U{(ji'~,) = {^l, ), with U being the negative of the matrix shown in Eq. ( 10.23 ), 



takes the form (01,) cx: (vi 
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is furthermore symmetric under the Klein symmetry generated by S and U of Eq. ( |10.22| ) . 
This can be shown either explicitly by calculating S^M^rS = U^MrjiU = Mrr^ or by 
realising that the alignments of the two neutrino flavons cf)-, and of Eq. ( |10.29| ) are 

3 3 

left invariant under both S and U (in the appropriate representations). With tulr being 
proportional to the identity matrix, the effective light neutrino mass matrix obtained 
from the type I seesaw formula, is symmetric under S and U as well. This particular Klein 
symmetry which originates from the family symmetry A(96) corresponds to the so-called 
bi-trimaximal mixing pattern [31,104] in the neutrino sector, 
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The bi-trimaximal mixing structure, which is a special form of TM2 mixing, does not 
involve any CP phases and translates to the following values of the neutrino mixing angles 
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sm 



-^(a_) « 12.2° , = 023 = tan"^ ( ^ - l) ^ 36.2° . (10.35) 



It is a remarkable fact that 9'(^ deviates from the experimentally measured value of the reac- 
tor angle ^13 ~ 9° by about 3°, a deviation which is typical for charged lepton corrections to 
a vanishing 13 neutrino mixing angle in SU{5) GUTs with implemented Georgi-Jarlskog 
mechanism. Indeed, taking into account the left-handed charged lepton 12 mixing with 



612 ~ A/3 and phase 612, see Subsection 3^, leads to the relation 

sin 013 w a_ - ^6)^2 cos 6I2 ■ (10.36) 

As the charged lepton correction to ^13 has to be negative in order to hit the measured 
value, one must assume 6^2 ^ 0. This entails a vanishing Dirac CP phase 

5^0 , (10.37) 

as well as the following phenomenologically viable lepton mixing angles, 

9i3 « 9.6° , 012 ~ 32.7° , 023 ~ 36.9° . (10.38) 

The charged lepton corrections are thus crucial for direct models based on a A (96) family 
symmetry. Furthermore, a specific choice of the phase 5^2 is required, which must eventually 
be explained in a more complete model, for example along the lines of the models proposed 



in [103], see also end of Subsection 7.2 



11. Conclusion 

Neutrino physics has progressed at a breathtaking rate over the last decade and a half, 
as discussed in Section ||, with a major discovery almost every other year, as can be seen 
by glancing at the milestones in Subsection |1.2| . The big experimental result of 2012 is 
the measurement by Daya Bay and RENO of the reactor angle, which has had a major 
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impact on neutrino physics, making the early measurement of CP violation possible. It 
has also ruled out a large number of models which predicted the reactor angle to be zero, 
including simple patterns of lepton mixing known as BM, TB or GR which are presented 
in Subsection |3.1| . As discussed in Section ^, these simple patterns can result from discrete 
symmetries, reviewed in Section |5|, such as A^, 84^ or A^. On the other hand, anarchy 
expected the reactor angle to be rather large and this is what was observed. Does this 
mean that we should give up on the symmetry approach in favour of anarchy? 

In this review the idea of giving up on the symmetry approach in favour of anar- 
chy is given short shrift for good reason, namely some simple mixing patterns remain a 
good approximate characterisation of lepton mixing. For example, we have seen that it 
is possible to describe the current pattern of lepton mixing by simply perturbing around 
TB mixing. This perturbation can be parameterised in terms of the reactor (r), solar (s) 
and atmospheric (a) deviations from TB mixing introduced in Subsection |3.2| , where such 
perturbations are no larger than the Wolfenstein parameter A, which apparently coincides 
with the reactor parameter r. Indeed TBC mixing provides a good approximation to the 
observed lepton mixing. 

From the symmetry perspective, the measurement of the reactor angle has caused 
theorists to work harder to explain the observed deviations from TB mixing, which in 
the near future could also include the atmospheric deviation parameter a and the solar 
deviation parameter s. Indeed there are already hints from the global fits to oscillation 
data that both these parameters could be non-zero, as discussed in Subsection |2.6| . The 
models based on small discrete family symmetries such as ^4, S4 or A^ may be viewed as 
predicting BM, TB or GR lepton mixing only at the LO in the absence of HO operator 
corrections or charged lepton corrections, not to mention renormalisation group running or 
canonical normalisation corrections. Indeed there are many sources of corrections that can 
modify the naive simple mixing patterns apparently predicted by discrete family symmetry. 

If we are very lucky then it is possible that only a special subset of these corrections 
are important, as discussed in Section ^. For example, if only Cabibbo-like charged lepton 
corrections are important then this leads to solar mixing sum rules. Although, in the 
framework of GUT models, it might seem natural to equate charged lepton corrections to 
the Cabibbo angle, this assumption is at odds with the simplest type of relations between 
charged lepton masses and down-type quark masses, where such relations would prefer 
the charged lepton mixing angle to be about a third of the Cabibbo angle. Alternatively, 
in the absence of charged lepton corrections, if only certain kinds of HO corrections are 
important then TB mixing could be reduced down to TMl or TM2 mixing, where only half 
of the original Klein symmetry of the neutrino mass matrix is broken, and this leads to 
atmospheric sum rules. In either case, solar and atmospheric sum rules provide interesting 
relations involving the deviation parameters r,s,a together with cos 5 which can be tested 
in future neutrino oscillation experiments. 

We have distinguished between two general approaches of using discrete family sym- 
metry to build realistic models, referred to as direct or indirect. In the direct approach 
the Klein symmetry of the neutrino mass matrix, discussed in Subsection is identified 
as a subgroup of the discrete family symmetry (with a different subgroup preserved in 
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the charged lepton sector), while in the indirect approach the Klein symmetry emerges 
accidentally after the discrete family symmetry is completely broken. Both direct and in- 
direct approaches rely on spontaneous breaking of the discrete family symmetry via new 
scalar fields which develop VEVs, referred to here as flavons, as discussed generally in 



Subsection 6.1 



A common feature of both the simplest direct or indirect models is the use of the type I 
seesaw mechanism, reviewed in Section with form dominance. In fact the different 
types of seesaw mechanism are generally reviewed along with the sequential dominance 
mechanism. As discussed in Subsection form dominance corresponds to the columns 
of the Dirac mass matrix being proportional to the columns of the PMNS matrix in the 
basis where the right-handed neutrino mass matrix is diagonal. The virtue is that this 
usually leads to a form diagonalisable physical neutrino mass matrix, with mixing angles 
being independent of neutrino masses. In the case of the simplest indirect models, form 
dominance is usually identified with constrained sequential dominance featuring a normal 
mass hierarchy with the lightest neutrino mass being approximately zero. The downside 
of form dominance is that, since the columns of the Dirac mass matrix are orthogonal, 
leptogenesis is identically zero so corrections to form dominance may be required. 

The direct approach is discussed in detail in Section |^, including flavon vacuum align- 
ment and model building strategies following Daya Bay and RENO. In the direct approach 
it is possible to achieve the simple patterns of lepton mixing, namely BM, TB or GR as 
an approximation to lepton mixing using small discrete family symmetries such as A4, S4 
or A^, and then consider the effect of corrections as discussed above. If we are lucky then 
such corrections could either preserve the Klein symmetry in the neutrino sector and break 
the symmetry in the charged lepton sector, leading to solar sum rules, or preserve half the 
original Klein symmetry in the neutrino sector, leading to TM mixing and atmospheric 
sum rules. Indeed it is possible to start with only half the original Klein symmetry in 
the neutrino sector arising as a subgroup of the family symmetry, as in for example, 
which we refer to as the semi-direct approach. More generally, however, the whole original 
Klein symmetry will be broken by HO corrections, and also charged lepton, renormalisa- 
tion group and canonical normalisation corrections will also be present. Alternatively it 
is possible to use larger discrete family symmetries in which the Klein symmetry already 
gives a non-zero reactor angle at the LO. For example A(96) can give bi-trimaximal mixing, 
where the reactor angle is non-zero and the solar and atmospheric angles start out equal, 
with all angles receiving modest charged lepton corrections. 

An analogous discussion for the indirect approach is provided in Section ^. An im- 
portant difference is that TB mixing at the LO can be achieved not only from groups 
in which the Klein symmetry Z2 x Z2 can be embedded but also in groups of odd order 
such as A(27) or Ty, or more generally infinite classes of groups. The types of correction 
to TB mixing discussed above are also possible for the indirect case, perhaps leading to 
solar and atmospheric sum rules if we are lucky. However, the indirect approach offers new 
alternatives to achieving a reactor angle already at the LO, without resorting to large dis- 
crete family symmetry groups, by using new vacuum alignments to construct the neutrino 
mass matrix. Such new alignments are more easily achieved with smaller groups since the 
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discrete family symmetry is completely broken and we are freed from the requirement of 
identifying the Klein symmetry as a subgroup. Examples of new vacuum alignments, which 
can be achieved even in A4, include PCSD and CSD2, where PCSD can give the successful 
TBC mixing if the misalignment parameter is identified with the Wolfenstein parameter. 

We have already explained that we prefer the symmetry approach to, say, anarchy, since 
simple symmetric mixing patterns still provide a good approximation to reality. However 
there is a deeper reason why we prefer to use symmetry, namely to address the flavour 
problem. In our view, to abandon the symmetry approach would completely miss the 
opportunity provided by neutrino physics of elucidating the flavour problem, the problem 
of all quark and lepton masses and mixing parameters, including CP violating phases. The 
history of physics, if it tells us anything at all, teaches us that symmetry and unification 
have always provided a guiding light in the path to understanding deep problems in physics. 
Therefore in this review we have considered models based not only on discrete family 
symmetry, but also using GUTs, in order to address the flavour problem. Motivated by 
such considerations, in Section ^ we have briefly review grand unified theories and how they 
may be combined with discrete family symmetry to describe all quark and lepton masses 
and mixing, tabulating some recent examples of this kind. 

Finally in Section ^ we have discussed three model examples which combine an SU (5) 
GUT with the discrete family symmetries A4, S4 and A (96). These models are presented 
in sufficient detail to illustrate the complexity of the current state of the art of GUTs of 
flavour that is required to account for all quark and lepton masses and mixing. Critics will 
use these examples as evidence that the effort of constructing such models is not worth the 
trouble and question what has been achieved by all this complexity. They will also point 
out that the number of input parameters in these models exceeds the number of mass and 
mixing parameters that we are trying to explain. However this last observation misses the 
point. What is relevant is the number of predictions (or postdictions) not the number of 
parameters. Typically each of these models contains half a dozen relationships such as the 
GST and GJ relations which agree with experiment, and neutrino mass and mixing sum 
rules giving predictions for future neutrino experiments. Thus one feels that something has 
been understood by constructing these models, including the mass hierarchy and origin of 
the three families as well as the milder neutrino hierarchy with bi-large lepton mixing. 

The flavour problem is not going away, in fact since 1998 it has got significantly worse 
with at least seven new parameters arising from the neutrino sector on top of the three 
charged lepton masses and the ten flavour parameters from the quark sector. However the 
neutrino parameters also provide some clues such as small neutrino masses, bi-large mixing 
and a Cabibbo-like reactor angle, and there is a ghost of a chance that these clues may be 
just enough to allow us to unlock the whole flavour puzzle. We are not there yet, but the 
hope is that the details of the admittedly rather complicated models given in this review 
may inform and inspire new young researchers to do better. 

It is still not too late for theorists to redeem their past failures to successfully predict 
anything in the neutrino sector by making a genuine prediction which can be tested by 
experiment. Indeed there is still room to make predictions for the solar and atmospheric 
deviation parameters s,a as well as cos 5, or to relate them via sum rules to the reactor 
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parameter r since all these parameters can and will be measured in future high precision 
neutrino experiments. Also one can make predictions for the pattern of neutrino masses 
including their ordering and their scale. A really crucial question is whether neutrinos are 
Dirac or Majorana particles, in the former case they would be just like the quarks and 
charged leptons, while in the latter case something qualitatively new is involved such as 
the see-saw mechanism for example. In the absence of experimental information about 
these questions, we must admit that we do not yet understand the origin of neutrino mass, 
which remains one of the biggest unsolved mysteries of the Standard Model. 
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A. Equivalence of different parametrisations 



In this appendix we exhibit the equivalence of different parametrisations of the lepton 
mixing matrix. A 3 x 3 unitary matrix may be parametrised by 3 angles and 6 phases. We 
shall find it convenient to parametrise a unitary matrix hy}^ 



V'^ — -P2-R23-Rl3-Pl-Rl2-P3; 



(A.l) 



where Rij are a sequence of real rotations corresponding to the Euler angles Qij , and Pi are 
diagonal phase matrices. The Euler matrices are given by 
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sin 9ij. The phase matrices are given by 
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By commuting the phase matrices to the left, it is not difficult to show that the 

(A.4) 



parametrisation in Eq. (A.l) is equivalent to 

W = PU23Ui3Ui2, 



where P = P1P2P3 and 
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^*It is convenient to define the parametrisation of rather than V because the lepton mixing matrix 
involves and the neutrino mixing angles will play a central role. 
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where 



S23 = X + ^2 - W3, 
^13 = wi - 003, 

612 = bOl — UJ2- 



(A.6) 



(A.7) 



(A.8) 
(A.9) 
(A.IO) 



The matrix U is an example of a unitary matrix, and as such it may be parametrised by 
either of the equivalent forms in Eqs. ( A.!] ) or (A. 4). If we use the form in Eq. ( A.4D then 
the phase matrix P on the left may always be removed by an additional charged lepton 
phase rotation AT4^ = P\ which is always possible since right-handed charged lepton 
phase rotations can always make the charged lepton masses real. Therefore U can always 
be parametrised by 

U = U23UUU12, (A.ll) 

which involves just three irremovable physical phases 5ij. In this parametrisation the Dirac 
phase 6 which enters the CP odd part of neutrino oscillation probabilities is given by 



<5 = (5i3 - 623 - <5l2- 
Another common parametrisation of the lepton mixing matrix is 

U = R23U13R12P0, 



(A.12) 



(A.13) 



where 

/ e^^i \ 

Po = e^^2 , (A. 14) 

V 1; 

and in Eq. ( [A.13 ) C/13 is of the form in Eq. ( [A.6D but with replaced by the Dirac 
phase 6. The parametrisation in Eq. (A.13) can be transformed into the parametrisation 
in Eq. ( |A.ll ) by commuting the phase matrix Pq in Eq. ( A.13D to the left, and then 
removing the phases on the left-hand side by charged lepton phase rotations. The two 
parametrisations are then related by the phase relations 

^^23 = /32, (A.15) 
Si3 = 6 + Pi, (A.16) 
<5i2 = /3i-/32. (A.17) 

The use of the parametrisation in Eq. ( [A.13 ) is widespread in the literature, however it is 
often more convenient to use the parametrisation in Eq. ( A.ll ) which is trivially related 
to Eq. ( [A.13 ) by the above phase relations. 
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B. Deviations from TB mixing to second order 
B.l PMNS matrix expansion 

The PMNS matrix when expanded to first order in the three real parameters r, s, a defined 
in Eq. 1^) reads [64] 



C^PMNS 



^(1 



VE 
\ 75 



A.{l + s-a + re^') ^(1 



73(1 + 

1 „ 



'1 + s + a- 



re 



V3' 
1 



ire 



75^""' 
73(1 + 



Pi 



12- 



+ a + ^(1 - a) J 



The second order corrections to the PMNS matrix elements are [64], 



(B.l) 



AUel 

AUe2 
AC/e3 

AUrl 
AUr2 
AUr3 



/2/ 1 2 ^ 2\ 

0, 

-^{-rse'^ 
-^{--rse^' 

J_(\rse'^ + rae'^ + so), 
v6 2 



-rae H — so 
2 2 



1 ,1 



^/3'2 



1(5 



-rae 



-sa 



-— r 
4 



The Jarlskog CP invariant to second order is then given by [64] 

,r rs. . ^ 

J PS (- H ) sm6. 

^6 12^ 

B.2 Neutrino oscillations in matter 



(B.2) 



(B.3) 



In this appendix we present the complete formulae for neutrino oscillations in the presence 
of matter of constant density to second order in the quantities r, s, a and A21, where it is 
assumed that A21 <^ 1. Here Ajj = 1.27AmfjL/E with L the oscillation length in km, E 



the beam energy in GeV, and Am~ 

VL 
2A 

V ^ 7.56 X lO^^V 



Am; 



A = where V is the potential expressed in units of eV as 



m?- in eV^. We write A = A31, a — . ^ 



^ and 



(B.4) 



where p is the matter density of the Earth in units of g/cm^ and A'^e ~ 0.5 is the number 
of electrons per nucleon in the Earth. The constant density approximation is good when 
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the neutrino beam only passes through the Earth's crust where p ~ 3 g/cm^ or the Earth's 
mantle where p ~ 4.5 g/cva? . 

The complete set of neutrino oscillation probabilities for electron neutrino or muon 
neutrino beams in the presence of matter of constant density to second order in the pa- 
rameters r, s,a and a are [64]: 

8 osin^^A osin2(.4-l)A , , 

4 osin^ylA nsin2M-l)A 
P — -rv \-r - 

^^"9 ^2 (A -1)2 

4 ^, sin AA sin(A - 1) A ,^ ^, 

+ -mcos(A - 6)^ P^^l)^' ^ ^ ^ 

4 ,sin2ylA ,sin2M-l)A 

'""9 (A -1)2 

4 ,. sin^lAsinfA- 1)A _ 
--r«cos(A-<5)-^^ (B.7) 

_4 2sin2ylA 2sin^(^-l)^ 

= 9« ^2 (A - 1)2 

4 ,,sinAAsin(A- 1)A 
+ -racos(A + ,5)— p^^' ^ 
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P^^ = 1 - (1 - 4a2) sin2 A + - (1 - s)aA sin 2A 

o 

4 2 sin2 AA 4 2^2 

a 7^5 q;2A2cos2A 

9 ^2 9 

4 2 1 / . sinAA , ^ , , A . 
+ qQ;2— smA — - — cos(yl-l)A- — sm2A 

2 sin2(A-l)A 

' (A -1)2 

1 .2 A ..... A sin(A - 1)A A 



sin A cos A A — A sin 2 A 



- 1 V - 1) 2 

4 ^ ^ sin AA sin(A - 1)A 
- -racos^cosA— {A-\) ' ^ 

2 4 
P^^ = (1 - 4a^) sin^ A - - (1 - s)aA sin 2A + -a^ A2 cos 2A 

o 9 

4,1/ , sin AA , ^ , ^ A \ 
- qQ:2— I sin A — - — cos(A - 1)A - — sin2A 

1 1 ( . ^ .A sin(A - 1)A A ^ . 
+ -r^ smAcosAA — j- -y — Asm2A 

4 . , sinylAsin(A - 1)A 
+ -rasm<5smA— Cl^^T) ' 
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C. Generators and Clebsch-Gordan coefficients of S4, A4 and T7 



In this section we list the generators of the groups 5*4, A4 and T7 in the basis where the 
order three generator is diagonal. As this basis is particularly convenient for model building 
purposes, we state the corresponding (basis dependent) Clebsch-Gordan coefficients for all 
non-trivial Kronecker products. We first consider the two intimately linked groups S4 and 
A4, before discussing the group Tj. 



C.l The groups S4 and A4 



The triplet matrix representations of the three S4 generators S, T and U in the T-diagonal 
basis can be obtained from Eq. (^.g|). Noticing that the h generator (corresponding to U) 
in Eq. (^^) occurs only quadratically, we immediately find another triplet representation 
by changing the overall sign of U . The obtained irreducible representations are called 3 
and 3', respectively. Likewise we find the two singlet representations 1 and 1'. Summing 
up the square of the dimensions of these representations, 1-^ -|- 1^ -|- 3^ -|- 3^ = 20, shows 
that there exists only one more irreducible representation, namely the doublet 2. Its 
matrix representation is presented, together with the other irreducible representations in 
the following table. 



1,1' 
2 

3,3' 



s 



u 




±1 




The same table also shows the representations of the 5*4 subgroup A4 , generated by S and 
T only. Dropping the U generator, it is clear that both triplets of 5*4 coincide with the 
single A4 triplet. Likewise, the two S4 singlets correspond to the trivial singlet of A4. The 
S4 doublet, on the other hand, becomes reducible once the U generator is removed. Hence, 
it decomposes into two separate non-trivial irreducible representations of A4, 1" and 1'. 

The non-trivial 5*4 product rules in the T-diagonal basis are listed below, where we use 
the number of primes within the expression 



a 



(') 



(3 



(') 



7 



(') 



(C.l) 
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to classify the results. We denote this number by n, e.g. in 3 3' ^ 3' we get n = 2. 



lW(g) 2 



iW^sW ^ 3(0 



1 (8)1 ^ 1 
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3' 
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3' 

3 
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3' 



02/32 + "3/51 + ai/53 
J-ir(a3/33 + ai/32 + a2^i), 

^2qi/3i - Q2/33 - 03/32^ 
203/33 - ai/32 - a2/3i | , 
^2a2/32 - "3/31 - "1^3, 

'^02/33 - a3/32 
ai/32 - a2/3i 
yas^i - ai/33^ 



The Clcbsch-Gordan coefficients can be obtained from these expressions by simply 
dropping all primes and identifying the components of the 5*4 doublet 2 as the 1" and 1' 
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representations of A4. We thus find the non-trivial ^4 products, exphcitly, 

1 

l'(g)3 3 

l"(g)3 ^ 3 

3 3 
3(^3 
303 



1 

1' 
1" 

3 + 3 



ai/3i + 02/33 + 03/32 
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Ol/32 
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03/32^ 

02/3l 

01/33/ 



C.2 The group T7 

The Frobenius group T7 = 
presentation, see e.g. [81], 



Zj XI Z3 is obtained from two generators a and c obeying the 



< a, c a 



1 , aca 



-1 



(C.2) 



Notice that with A: = 2 in Eq. ( |5.11D , the condition 1 + 2 + 2^ = mod 7 holds. A 
triplet representation with non-diagonal order three generator is given in Eq. ( |5.12| ), where 
r] = e~ . To diagonalise a, we apply the basis transformation of Eq. ( |5.10| ) (followed by 
the phase transformation in order to bring the c generator into a more appealing form), 
resulting in the matrix representation of the triplet 3 as shown in the below table. Clearly, 
also contains another triplet representation given by the complex conjugate 3 of the 3. 
Furthermore there are three singlet representations. 
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Although the order-seven generators of the triplet representations look rather involved, the 
Clebsch-Gordan coefficients take a relatively simple form. Omitting the trivial products, 
i.e. those involving the singlet 1 as well as products of only one-dimensional irreducible 
representations, the product rules are reported below. Again, we use the convention that 
the components of the first representation of any given product a (8) /3 are denoted by 
while we use Pi for the components of the second representation. 



1' 



3 
3 

3( 
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